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ABSTRACT 

A computer simulation model is developed for use in 
analysing the Navy's Fleet Air Defense problem. The model 
provides for employment of monostatic and bistatic search 
radars in a clear or jamming environment. Use of the model 
is demonstrated by analysing a hypothetical problem. 
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I. INTRODUCTION 



The purpose of this thesis is to develop a computer 
simulation model for use in analysing the Navy's Fleet Air 
Defense problem, and to demonstrate the usefulness of the 
model with a hypothetical air defense scenario. The model 
is unique in that it provides for employment of bistatic as 
well as monostatic search radars. A bistatic radar is sim- 
ply a radar system in which the radar transmitter and the 
radar receiver are not collocated. This capability was pro- 
vided to allow for the investigation of the possible advan- 
tages of the use of a bistatic radar receiver in view of the 
extensive capability for passive detection of active radars. 

The model, acronym FADS, for Fleet Air Defense Simula- 
tion,' is programmed in FORTRAN IV level G, and is designed 
to run on an IBM 360/65 computer meeting the core require- 
ment of 350 kilobytes. 

FADS is structured along the lines of FLOATS IIB, a 
model currently maintained at the Johns Hopkins University 
Applied Physics Laboratory [Ref. 6]. The significant dif- 
ferences from FLOATS IIB which have been incorporated in 
FADS are: bistatic search radar systems, provision for en- 

gagement of bombers by surface-to-air missiles (SAM), the 
addition of a surface-to-surface missile (SSM) capability 
for the defensive forces, a control logic which requires 
the defensive forces to hold fire until fii’ed on, the in- 
clusion of a dispersed defensive force disposition, and the 
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means for constructing up to six distinct target lists for 
the attacking units. 

FADS is suitable for analysing both large scale problems 
and one-on-one situations in a clear or jamming environment. 
In general, the model provides for the description of up to 
90 enemy raiders of 10 different classes, 40 defensive units 
of 10 different classes, 15 classes of radars, 10 classes of 
SAMs, 35 time delay distribution functions, and 18 logic 
switches are included which vary the defensive play of the 
simulation. During the computer debugging phase of the 
model construction, FADS was exercised in all of its various 
modes. Included among these were clear and jamming environ- 
ment, single and multi threats, area and point defense, 
aim monostatic ano. oxstatic search radar systems. 

The particular problem used to demonstrate the useful- 
ness of the model is an investigation of the relative ef- 
fectiveness of 100 nautical mile and 50 nautical mile 
non-fire control dedicated SAM systems deployed on picket 
type surface units with bistatic search radar capability. 
This type of SAM system employs active or passive homing 
throughout its entire flight. The measure of effectiveness 
for the respective SAM systems is the rate at which enemy 
air-to-surf ace and surface-to-surface missiles (ASM/SSM) 
arrive at their target area during the ASM/SSM impact phase. 
This is an indication of the level of the enemy attack sur- 
viving the picket defenses and will demonstrate the dimin- 
ishing of the intensity of the air battle which will be 
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waged in the vicinity of the main force. The primary concern 
is to determine what effect the different missile ranges 
will have on breaking up this focus of an enemy ASM/SSM 
attack . 

The analysis was conducted with hypothetical or illus- 
trative values for all systems. Specific performance char- 
acteristics for the bistatic radar system and the missile 
system do not exist. Care was 'taken to ensure that the hy- 
pothetical values were within the range of existing systems 
and that the proposed systems' parameters were comparable to 
present system capabilities when appropriate. 

Section II of the thesis contains a general description 
of the flow of events in FADS and a discussion of the princi- 
ple assumptions in the model. The specific inputs and re- 
sults of the demonstration problem together with an analysis 
of the results are contained in Section III. Sections IV 
and V provide a detailed listing of the inputs required and 
the outputs received from FADS. Conclusions are given in 
Section VI. The appendixes contain background information 
on bistatic radar systems, a verbal description of the simu- 
lation program, glossaries , and a FORTRAN IV computer program 
listing. 
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II. BACKGROUND AND DISCUSSION 



A. GENERAL DESCRIPTION OF THE SIMULATION MODEL 

FADS is a large scale store-event simulation of a Task 
Force size surface missile engagement with an enemy air 
attack. The model is capable of representing a threat con- 
sisting of a coordinated missile attack launched from air- 
craft, surface units or submarines. Although the model 
provides for defensive engagement of the enemy missiles and 
launch platforms with either SAM or SSM, there is no deter- 
mination of defensive interceptor attrition of enemy systems 
during the engagement. If interceptors are to be employed, 
they must be played against the launch platforms prior to 
the start of the game and the enemy forces reduced to the 
appropriate level. The Air Battle Analyzer developed by the 
Applied Physics Laboratory of Johns Hopkins University is a 
suitable tool to determine these interceptor attritions 
[Ref. 11] . 

The simulation determines friendly force equipment avail 
abilities at the start of each iteration. These availabil- 
ities are based on a comparison between user inputs of the 
probability that a given equipment is operational and the 
value of a uniform (0,1) random number which is drawn by the 
computer. The next step in the simulation's iterative proc- 
ess is to determine detections and firing events for the 
friendly forces against the enemy launch platforms. These 
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SAM/SSM fire events are stored in an events list in a time 
sequence. The enemy ASM/SSM launch sequence and targeting 
is computed next. Detections and resultant fire events 
against the enemy ASM/SSM' s are then determined and entered 
in the appropriate sequence in the events list. Friendly 
SAM/SSM intercepts and target hits are evaluated followed by 
enemy ASM/SSM hits and ship kill determinations. The kill 
and damage events are time interactive so that offensive and 
defensive force capabilities are decremented as hits occur. 
The simulation then prints statistics for the iteration and 
starts the next iteration by re-entering the flow at the 
point where equipment availabilities are determined. The 
program stops upon completion of the specified number of 



iterations and cumulative statistics arc provided. 



B. ASSUMPTIONS 

1. Bistatic Radar System 

In the bistatic detection portion of the simulation, 
the principle assumptions made regarding bistatic radar sys- 
tems are: a) the target's bistatic radar cross section is 
equal to its monostatic cross section; b) the target is not 
located along the transmitter-receiver base line; and c) the 
radar propagation losses are the same as those experienced 
in monostatic systems. 

a. Radar Cross Section 

The approximation for the bistatic radar cross 
section is made on the basis that, in the limit, the bi- 
static cross section is equal to the rnonostatic cross 



15 



section except for the case where the scattering angle is 
equal to 180° [Ref. 10, p. 590]. The scattering angle is 
the angle between the path from the transmitter to the tar- 
get and the path from the target to the receiver. It is 
equal to 180° when the target is on the base line between 
the transmitter and the receiver. The fact that the mono- 
static and bistatic cross sections are the same in the limit 
does not imply that they are everywhere equal. It does 
imply that the range of the values will be equal and thus on 
the average allow usage of the monostatic cross section to 
provide representative values for detection ranges during 
the iterative process of the computer simulation. 

b. Target Location 

The assumption concerning the target not being 
on the base line between the transmitter and the receiver is 
necessary since range ambiguities result when the scattering 
angle is equal to 180°. This assumption is not unreasonable 
since the typical track may cross the base line, but will 
not originate on it and remain on the base line. A situa- 
tion in which this could occur would be the ca.se when a 
raider while on the base line launches a missile at a target 
which is either the transmitter or the receiver platform. 

c. Propagation Losses 

The assumption concerning the propagation losses 
is appropriate in that similar assumptions are routinely 
made in arriving at detection ranges for monostatic radars. 

A representative propagation loss for the nominal detection 
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range is used in conjunction with the signal to noise ratio 
required for a specified probability of detection to arrive 
at a minimum reflected power level which must be received 
in order for target detection to occur. 

An additional assumption concerning the bistatic 
radar system in the simulation model is that the bistatic 
transmitter is always in an operative status. It is not de- 
stroyed by enemy missiles during .the attack phase nor is it 
evaluated as being down or degraded during the equipment 
status evaluation which is conducted at the start of each 
iteration. Individual bistatic receivers have an assigned 
probability of being operational which should also reflect 
the reliability of their associated transmitter. 

2 . Jamming 

Some major assumptions were made in simulating the 
jamming environment in order to keep the problem tractable 
and to remain within reasonable time limits for the computer 
program. The jamming assumptions will tend to drive the 
solution to conservative estimates for the defensive forces. 
The first assumption made is that the jammer and detecting 
radar have a single gain for their main beam. Thus targets 
not in the main beam are considered to be in a unique side- 
lobe which is a specified number of db below the main beam. 
To determine the level of the jamming power arriving at the 
detecting radar, it was necessary to approximate the posi- 
tion of the jammer relative to the detecting radar. This 
was done by determining the jammer's position at the time a 
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given target was at the midpoint of its track through the 
radar's detection envelope. In the case of monostatic radar 
this point was the target's closest point of approach. For 
bistatic radar the midpoint in the first detection envelope 
was utilized. The jamming power arriving at the detecting 
radar, given the geometry of this midpoint position, was 
then used to calculate the detection range for the target. 
Additionally, in the computation of bomber detections and 
"shoot" events, no provision is made for jammers which may 
have been shot down prior to the time of detection. This 
assumption is not required in the case of missile detections 
where detection events for missiles take into consideration 
reduction in jamming levels caused by previous destruction 
of jammers. 

The model is limited to four frequency bands for 
frequency separation of individual radars in the jamming en- 
vironment . 

3 . Target Tracks . 

The principle target track assumptions are concerned 
with minor anomalies from the real world situations. These 
assumptions do not detract excessively from the simulation 
results, but do ease computation and simplify programming. 

A fixed altitude for the raider and missiles, and the lack 
of course changes on the part of raiders are such simplify- 
ing assumptions. In the case where a raider or missile 
track does not change more than 0.1 nautical mile (nm) in 
either the x or y direction, the target is assumed to fly 
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parallel to the non-changing axis. If the target does not 
move more than 0.1 nm in both the x and y directions, the 
target is assumed to be stationary. Additionally, raiders 
who have completed their missile launches continue on their 
preset course rather than exiting the battle area. Although 
these targets will be engaged by the defensive forces if 
higher priority targets are not present, the resulting tar- 
get hits can be readily "backed out" of the simulation re- 
sults by the user if desired. On the other hand, this 
assumption does allow raiders that are providing standoff 
jamming protection for the missiles to continue to close the 
defensive forces and thereby increase the effectiveness of 
their jamming. 

Ocher less general, assump cions are covered as they 
occur in the details of the subroutine descriptions in Ap- 
pendix B. 
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III. 



DEMONSTRATION PROBLEM 



This section presents a demonstration of the use of FADS 
in analysing a hypothetical problem involving bistatic radar 
systems and surface-to-air missiles (SAM). The problem is to 
select the maximum range required for the SAM system to be 
employed by a small surface combatant (SSC). The proposed 
SAM system is to protect the major units of a naval task 
force involved in an engagement with enemy missile systems. 
The SSC will have a bistatic search radar and the SAMs will 
be active homing missiles requiring no fire control director. 
As such, the SSC will emit no radiation. The SAM character- 
istics and performance capabilities have been specified with 
the exception of the maximum range. Two range capabilities 
will be considered, a 100 nautical mile missile and a 50 
nautical mile missile. The problem thus becomes that of de- 
termining the relative effectiveness of the SSC when equipped 
with the different range missiles. As discussed in the in- 
troduction, the measure of effectiveness used to evaluate 
the missiles is the number of enemy missiles penetrating to 
the force center per unit of time. That is, the number of 
missiles arriving at the force center per minute during the 
time period in which missiles are impacting at the force 
center. 



20 



A. THE SCENARIO 



The BLUE force consists of one High Value Target (HVT) 
and eight small escorts. The BLUE force is maneuvering 
within range of enemy air forces. The HVT is located at 
force center and the pickets are equally spaced on a circle 
one hundred miles from the HVT. The picket force is silent 
except for one early warning aircraft which carries the 
transmitter for the bistatic radars. For purposes of this 
analysis, the aircraft is assumed to be at force center and 
at 20,000 feet. 

The RED force has detected the radar transmissions of 
the aircraft and the HVT and have launched an attack of 
twenty-one aircraft accompanied by three jammers. The RED 
aircraft will attack siinui Laneously in three sections, 
separated by forty degrees of bearing. The RED air force 
has practiced this maneuver to the point that the maximum 
deviation from simultaneous launch will be fifteen minutes. 
Each aircraft will take a maximum of three minutes to launch 
its missiles. The three sections will approach their target 
from 320, 000, and 040 degrees true. The jammers will carry 
equipment to cover all radar bands used by the BLUE force. 
Each section will consist of one heavy and six light attack 
aircraft . 

The RED air force will attempt to penetrate the BLUE de- 
fense and sink the HVT. The HVT is considered to be unarmed. 
In order to target a BLUE unit in this dispersed formation, 
a RED unit will have to have positive identification of the 
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unit as hostile. An ECM fix, or visual ID will be adequate. 
As the BLUE pickets are not radiating, no ECM fix may be ob- 
tained. As the pickets are assumed to be smaller than a 
destroyer escort, and the RED aircraft will fly at a minimum 
altitude of 33,000 feet, no positive visual identification 
will be allowed. 

The radio and radar silence of the pickets leaves them 
without the ability to coordinate their defense. No BLUE 
intercept will be allowed if the intercept point is within 
two miles of another BLUE unit's assigned position. The 
BLUE units will use a shoot- look-shoot doctrine. 

B. THE INPUTS 

1 . The Units 

a. The BLUE Units 

The BLUE units consist of the HVT, eight pickets, 
and one early warning aircraft. Positions are as indicated 
in the scenario. 

(1) HVT . For the purpose of this analysis the 
HVT is assumed to be unarmed. 

(2) Pickets . The pickets will employ a type 4, 
bistatic, band 4 radar with an antenna height of 50 feet. 

They will carry thirty missiles in canisters. The probabil- 
ity of a correct target identification will be .95. 

(3) Early Warning Aircraft . The BLUE early 
warning aircraft will be unarmed. It will carry a type 4 
radar and serve as the transmitter for the pickets' bistatic 
receiver . 
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b. The RED Units 

The three classes of RED units are heavy attack, 
light attack and jammers. 

(1) Heavy Attack . The RED heavy attack air- 
craft will carry a type 5 radar and have the following char- 
acteristics . 

WEAPONS: 1 type I ASM 

4 type II ASM 

SPEED: 415 Knots 

JAMMERS: none 

RADAR CROSS SECTION: 40 square meters 

ALTITUDE: 35,000 feet 

(2) Light Attack . The RED light attack air- 
craft will carry a type 5 radar and have the following char- 
acteristics . 

WEAPONS: 2 type II ASM 

SPEED: 420 Knots 

JAMMERS : none 

RADAR CROSS SECTION: 13 square meters 

ALTITUDE: 33,000 feet 

(3) Jammers . The RED jammer will carry no 
radar. It will have the following characteristics. 

WEAPONS : none 

SPEED:' 417 Knots 

JAMMERS: bands 1-4 

250 watts/MHZ /band 
359° beam width 
-20 db side-lobe ratio 
jamming target, HVT 
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RADAR CROSS SECTION: 32 square meters 

ALTITUDE: 40,000 feet 



2. The Radars 

The two radars which the simulation will use to play 
the simulation are the class 4 and class 5 radars. Class 4 
radars are found on all pickets and class 5 on all raiders. 

a. Class 4 Radar 

TYPE: Bistatic 

.1 

3: 120 nm/m 2 

x 

a: 30 nm/m 2 

SIDE LOBE RATIO: -20 db 

MAXIMUM INSTRUMENTED RANGE: 200 miles 

RADAR FREQUENCY BAND: 4 

BEAM WIDTH: 10 degrees 

TRANSMITTER UNIT: aircraft overhead HVT 

TRANSMITTER ANTENNA HEIGHT: 20,000 feet 

b. Class 5 Radar 

TYPE: Monostatic 

1 

3: 100 nm/m 2 

x 

a: 50 nm/m 2 

SIDE LOBE RATIO: -20 db 

MAXIMUM INSTRUMENTED RANGE: 200 miles 

RADAR FREQUENCY BAND: 3 

BEAM WIDTH: 3 degrees 
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CUMULATIVE PROBABILITY CUMULATIVE PROBABILITY 



3. The Reaction Time Distributions 



The reaction time distributions for target process- 
ing by the picket are displayed below. 




DETECTABILITY TO DETECTION DELAY IN MINUTES 

Figure 1 




DETECTION TO ID AND DESIGNATION TO WEAPONS EQUIPMENT 

Figure 2 
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CUMULATIVE PROBABILITY CUMULATIVE PROBABILITY 




ASSIGNMENT TO FIRING CONSOLE IN MINUTES 

Figure 3 




26 



4 . The Weapons 



a. The BLUE Weapons 

There are tv/o types of BLUE weapons. Both are 
SAMs and are denoted as SAM LR (long range) and SAM SR 
(short range). 

(1) SAM LR . This will be the 100 mile version 
of the BLUE missile. It will have two sets of intercept 
parameters. The first is for a large, slow, high flying tar- 
get. The second is for a small, fast, medium altitude target. 




SAM Intercept Contour for High Altitude Engagement 

The input value for the SAM LR class are: 
A = 100 
B = 100 
N = 2.0 
SLOPE = -20 
SEEK = -5 
SPEED = 1400 Knots 
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PROBABILITY KILL FOR SINGLE SHOT (PKSS) 

= .88 

PKSS (CROSSING) = .79 
SYSTEM RELIABILITY - .95 
LAUNCH RELIABILITY = .99 
MAXIMUM ALTITUDE = 59,000 feet 
MINIMUM RANGE =0.5 miles 
The intercept parameters A, B, N, SLOPE, and 
SEEK are defined in the description of subroutine TIMMY in 
Appendix B. 




Figure 6 

SAM Intercept Contour for the Medium Altitude Engagement 

A = 100 
B = 90 
N = 2.8 
SLOPE = -1 

All other parameters are the same as 
high altitude engagement 
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(2) SAM SR . This will be the 50 mile version 



of the BLUE missile. The intercept contours will be the 
same as those for the SAM LR except the down range and cross 
range values are half the SAM LR values, 
b. The RED Weapons 

The RED aircraft carry two types of ASMs. 

(1) ASM Type I . ' 

SPEED: 750 Knots 

ALTITUDE: 20,000 feet 

LAUNCH RANGE: 120 miles 

LAUNCH RELIABILITY: 0.85 

INFLIGHT RELIABILITY: 0.80 

RADAR CROSS SECTION: 1.79 square meters 

(2) ASM Type II . 

SPEED: 900 Knots 

ALTITUDE: 15,000 feet 

LAUNCH RANGE: 110 miles 

LAUNCH RELIABILITY: 0.85 

INFLIGHT RELIABILITY: 0.9 

RADAR CROSS SECTION: 1.01 square meters 
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C. RESULTS 



A summary of the results for 100 iterations is presented 
below. 



ASM Results. 


SAM SR 


SAM LR 


Available 


5100 


5100 


Destroyed with launch platform 


172 ( 3.4%) 


246 ( 4.8%) 


Launched 


4928 (96.6%) 


4854 (95.2%) 


Failed at launch 


724 (14.7%) 


750 (15.4%) 


Failed in flight 


429 ( 8.7%) 


418 ( 8.6%) 


Destroyed by SAMs 


0 


0 


Impact on HVT 


3775 (76.6%) 


3686 (75.9%) 


SAM Results. 


Available 


24000 


24000 


Launched 


3401 


6440 


Kills 


2116 (62.2%) 


2398 (37.2%) 


ASM 


0 


0 


Raiders 


2116 


2398 


Overkills 


605 (17.8%) 


2694 (41.8%) 


Failures 


680 (20.0%) 


1348 (20.9%) 


ASM Impact Intervals (minutes) . 


Mean 


1.468 


1.432 


Variance 


0.033 


0.037 


Skewness 


-0.060 


0.070 


Kurtosis 


-0.620 


-0.532 


Minimum 


1.091 


0.966 


Maximum 


1.849 


1.862 


ASM Impacts per Minute. 


Mean 


26.101 


26.192 


Variance 


14.651 


17.541 


Skewness 


-0.584 


0.453 


Kurtosis 


-0.006 


-0.136 


Minimum 


18.422 


18.129 


Maximum 


37.567 


37.246 



5. Histograms of the ASM impacts per minute are shown in Figures 7 
and 8. 
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Histogram of the ASM Impacts per Minute, SAM LR 

Figure 8 • 
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D. ANALYSIS 



The analysis of the results consists of three parts. 

The first part contains a parametric evaluation of the ASM 
impacts per minute for the two SAM systems. The second 
part provides a nonparametr ic evaluation. Part three pro- 
vides general comments concerning anomalies in the results. 

1 . Parametric Evaluation 

The purpose of the parametric evaluation is to test 
whether the mean ASM impacts per minute against the SAM SR 
(u 50 ) is significantly greater than the mean ASM impacts per 
minute against the SAM LR (u 100 ). With its shorter range 
capability, the SAM SR system should attain a lower number 
of kills than the SAM LR system. Consequently the number of 
ASM impacts on the HVT that occur during the defense by the 
SAM SR system should be greater than the number of impacts 
during the SAM LR defense. 

a. Null Hypothesis 

H 0 : u 50 £ u 100 ; SAM SR effectiveness is equal to 
or greater than the SAM LR system's as indicated by the mean 
ASM impacts per minute. 

b. Alternative Hypothesis 

II i : u 50 > u 100 ; SAM SR is less effective than 
SAM LR and therefore has a higher mean ASM impacts per min-, 
ute. 



H 0 . 



c. Significance Level 

A significance level of 0.05 was used in testing 
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d. Statistical Test 

The null hypothesis was tested against the al- 
ternative hypothesis using a t test. The t statistic for 
testing equal means is: 

x 5 o - x 1 0 o 

t = 

4 2 . 2 i 

r n 5 o +n i o o -I ■- n ^ o s n + n i o o s i 0 0 . 

l n 5 on 10 o 1 1 n 5 o + n] 0 o-2 J 

where : 



x 5 0 


= sample 


mean for 


SAM 


SR. 


(26.10) 


X 1 o 0 


= sample 


mean for 


SAM 


LR. 


(26.19) 


s ^ 

b 50 


= sample 


variance 


for 


SAM 


SR. (14.651) 


S 1 0 0 


= sample 


variance 


for 


SAM 


LR. (17.541) 



n 5 o“ n 10 o = number of samples for SAM SR and SAM LR. (100) 

The null hypothesis will be rejected if the computed test 
statistic t is greater than or equal to the table value for 
the 0.95 point of t distribution with 198 degrees of freedom. 
The table value is: 

t =1.96 

. 95,198 

e. Decision 

The computed test statistic is t = -0.1578. 

Since the computed t ratio is less than the tabulated value, 
the null hypothesis cannot be rejected at the 0.05 level. 

f. Conclusion 

The effectiveness of the SAM SR system was equal 
to or greater than the effectiveness of the SAM LR system at 
the 0.05 level of significance. 
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g. Justification fox' Use of the t Test 

In order to use the t test for difference in 
means, three basic conditions • had to be met. These were: 
i) the sample means were independent; ii) the underlying 
distributions from which the samples were drawn was normal; 
and iii) the underlying distributions had equal variances. 

(1) Independence . The independence of the 
sample means was a reasonable assumption since the individ- 
ual samples were obtained from different simulation runs 
using different random number streams. 

(2) Normality . The normality of the sample 
mean values could be asserted on the basis of the sample 
size (100) and the central limit theorem or by demonstrating 
that the means were computed from samples which came from a 
normal distribution. The underlying normal distribution can 
be shown on the basis that both samples were unimodal and 
symetric, and that the plots of their order statistics 
against the inverse normal pi'obability scale was linear. 

The unimodality of the samples is seen in the histograms in 
figures 7 and 8. Symmetry is evident in both histograms and 
supported by the fact that the kurtosis of the sample values 
is close to zero. The order statistic plots, which are 
clearly linear, are shown in figures 9 and 10. 

(3) Equal Variance . The relationship between 
the vai'iance of the population fi'om which the samples were 
drawn is tested using the F ratio. This test is applicable 
since the populations are normally distributed. 
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Inverse Normal vs SAM SR Order Statistics 
Figure 9 




Figure 10 
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(a) Null Hypothesis. H q : o* q = o| o(i ; the 
population variance (a| ) of the SAM SR system is equal to 
the population variance (a^ oo ) of the SAM LR system. 

(b) Alternative Hypothesis. Hj : o^ 0 f a 2 0 0 
the population variances are not equal. 

(c) Significance Level. A significance 
level of 0.05 was used in testing H Q . 

(d) Statistical Test. H Q was tested 
against IIj using the ratio of the unbiased estimators of the 
two population variances. 

F = —12. = 17 -' - — - = 1.1973 
S50 14 - 651 

The null hypothesis will be rejected if the computed ration 
F is greater than the tabulated value of F,, , _ „ . 

( i -a/2 ; 9 9 , 9 3 ) 



or less than F 


(a/2 ; 


9 9 , 


9 9)' 


The tabulated values are 


F < 


. 9 7 5 , 


99 , 


9 9) 


1.451 


F ( 


. 0 2 5 , 


9 9 , 


9 9) 


0.684 



(e) Decision. Since the computed ratio F 
is not greater than or less than the tabulated values for 
the F distribution, the null hypothesis cannot be rejected. 

(f) Conclusion. The samples were drawn 
from populations which had equal variances. 

2 . Nonparametr ic Evaluation 

The purpose of the nonparametric evaluation is to 
test whether the sample values for the SAM SR system and 
the SAM LR system are from the same population. 
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a. Null Hypothesis 



II 0 : There is no difference in the number of ASM 
impacts per minute when defending with the SAM LR system or 
with the SAM SR system. 

b. Alternative Hypothesis 

Hi : There is a difference. 

c. Significance Level 

A significance level of 0.05 was used in testing 
Ho . 

d. Statistical Test 

The Kolmogorov-Smirnov (K-S) two-sample test 
v/as used to test the null hypothesis. The K-S test requires 
the use of independent samples. Independence was shown in 



I 1 /IN b ▼ T » 

' • jl • £3 • ^ y 



tistic 



f n \ -y o T/' Q 4- -4- 4 ^ O 4- ^ #-1 ^ < r** 4 ^ 4" ^ ^ 4- ^ 

xilc uoou io uvjuuuo Leu unc bud” 



K 



max 



F (X) - G (Y) 
n v • m ' ' 



where F n (X) is the empirical cumulative distribution func- 
tion for the first set of sample values and G m (Y) is the 
empirical cumulative distribution function for the second. 

The null hypothesis will be rejected if there is 
a large enough deviation between the two sample cumulative 
distributions. The critical values of D for the level of 
significance of 0.05 and 100 observations in each sample is 
0.19233. 
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Decision 



e . 

The computed value of K was determined to be 
0.0800. Since the computed value of K is less than the 
tabulated critical value, the null hypothesis cannot be re- 
jected at the 0.05 level. A plot of the empirical cumula-- 
tive distribution functions for the two sets of sample 
values is provided in figure 11. This plot clearly shows 
the close relationship between the two underlying popula- 
tions . 

f. .Conclusion 

There was no difference in the defensive capa- 
bilities of systems in the environment in which they were 
examined. 

3 . General Comments 

The principal anomalies present in the data include 
the following: 

a) a lack of difference in the SAM LR and SAM SR 
effectiveness ; 

b) no SAM kills of ASMs; 

c) a greater number of overkills for the SAM LR; 

d) the number of SAM failures; and 

e) the short ASM impact interval. 

These anomalies do not significantly detract from the con- 
clusions presented in the parametric and nonparametric eval- 
uations . 
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Plot of Empirical Cumulative Distribution Function for Both SAMs 

Figure 11 



a. Lack of Difference in Effectiveness 

The most probable reason for the lack of differ- 
ence in the effectiveness of the two SAM systems was the 
jamming environment. The reduced detection ranges experi- 
enced in jamming will tend to preclude the use of SAMs at 
the longer ranges available to the SAM LR system. 

b. No SAM Kills of ASMs 



The lack of SAM kills of ASMs was primarily due 
to the geometry of the ASM launch positions and the program 
logic concerning SAM targeting. The program logic requires 
that the closest target be considered as most threatening 
and therefore engaged first. In the problem, the ASM launch 
positions were less than 10 miles from their point of clos- 



appro a C lx 



CAM ^ U -J ^ 



per 



• i rvn n 



Consequent ly , the 



ASMs with their greater speed were closer than the raiders 
to the SAM ships for only the first 30 seconds after their 
launch. After that period, their speed caused them to be at 
a greater range from the SAM ship. Thus the SAM ships had a 
30 second period in which to detect, acquire, and designate 
an ASM as the most threatening target. The reaction time 
distributions for the SAM ships precluded such a short re- 
action time in most cases, 

c. Overkills 

The higher number of overkills for the 100 mile 
SAM was due to the longer flight time of the SAM. This 
coupled with a lack of coordination resulted in a greater 
opportunity for another unit to obtain a kill on the target 
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while a second SAM was in flight. Targets that had been 
destroyed prior to the scheduled SAM launch time were not 
engaged. 

d. SAM Failures 

The number of SAM failures includes SAMs which 
were fired on the ASMs that impacted on the IIVT while the 
SAM was in flight. These attempted engagements occurred 
since the SAM units were not given the capability to deter- 
mine the ASM impact point. 

e. Short ASM Impact Interval 

The tightness of the ASM impact interval may 
appear to be pessimistic in relation to a real world threat. 
However, this should not detract from the comparison of the 

J QAU „ ^ J — 

IWU O tA ill £>y u . 



y 
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IV. MODEL INPUTS 



The inputs required to use FADS fall into three groups. 
These are the X vector, the target list assignments, and the 
target lists themselves. The amount of actual input re- 
quired depends upon the size of the problem to be analyzed. 

A. THE X- VECTOR 

All X vector inputs are real values and are read into 
the X vector in the following way. Each input card must 
contain six numbers. The first number is the integer address 
in the X vector. This integer should be placed in columns 
1-10, right justified. The remaining five numbers are the 
values to be inserted in the five X vector positions begin- 
ning with the position indicated in columns 1-10 of the input 
card. The format for each card with X vector data is 
(110, 5F10. 3) . If there are no data to be read into a posi- 
tion, no input is required, unless the position falls within 
the scope of another input card. The end of the X vector 
inputs is signaled by the inclusion of a zero card. This 
card contains a zero in column '10 and zeros in the five data 
positions. A brief overview of the organization of the in- 
put array is provided below. Details concerning the indi- 
vidual sub-sections are provided at the beginning of each 
sub-section. Specifics dealing with the input values are 
discussed in the description of the particular input value 
and as footnotes to the individual sections. 
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The following is the organization of the input X vector: 

X(l) to X(100), Control switches and program constants 

X(101) to X(500), Ship Classes 

X( 501 ) to X( 900) , SAM Classes 

X(901) to X(1200), Radar Classes 

X(1201) to X(1400), Ship Positions 

X(1401) to X(1800), Delay Functions 

X(1801) to X(2000), Raid Classes 

X(2001) to X(2200), Raid Weapon Classes 

X(2201) to X(2740), Initial Rajder Positions. 

1 . Control Switches 

These inputs include logic switches for varying the 
general play of the game. They also provide the means for 
constrolling the output in the form of debug and array 
printout switches. Several of the control switches are re- 
lated to inputs for specific systems which appear in later 
sections and must be consistent with later usage. Other 
portions of this input section are used to store values for 
generating statistics and for internal use by the model. 

These portions should not have inputs since their value will 
be set as required by the model. 

X( ) Inputs 

1 Case number 

2 Blue force equipment availability switch: 

0 = 100% availability 

1 = use input availabilities 
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M—2 



Inputs 



3 



4 



5 



6 



7 



8 



9 



10 



11 



12 



13 



Blue force radar detection range printout in 
Event 2: 

0 = no print 

1 = print 

Fleet coordination of Blue forces: 

0 = none 

1 = coordinated defense 

Blue ship out of action after N hits: 

0 = yes 

1 = no 

(N is a ship class input) 

Number of iterations. (Must be at least 1. If 
X(85) used, must be at least 10 and less than 
101 . ) 

Horizon/Sector Scan for Blue fire control radar 
system enabled: 

0 = no 

1 = yes if equipment capable 
Blue SAM Sector enable: 

D : — non ono-tn-po 

^ Uvu i/ Uo o »_> w u w a v_. 

1 = .sectors used 

Print control: 

0 = none 

1 = selected output tables and intermediate 

debug printouts. Set to 0 by the model 
after first iteration. 

Detectability to detection delay: 

0 = use input distribution function 

O.XX = a fraction; time for Raider or ASM to 
move this portion of calculated de- 
tection range. 

ISEED value. Initial seed value used in genera- 
tion of random numbers. Can be up to a seven 
digit number. 

X-array print out switch: 

0 = no print 

1 = print out X-array as read in, target 

lists, and target list assignments. 

Random number counter. Set to 0 by model. 
Maintains running count of random numbers drawn. 
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U 1 


Inputs 


14 


Total number of random numbers used. Computed 
by model. 


15 


Converted X- array print out switch: 

0 = no print 

1 = print out X- array as converted by model. 


16 


Blue first shoot control: 

0 = Blue may fire as soon as Red detected. 

1 = Blue may only fire after first Red 

missile launch. 


17 


Number of jamming raiders: 

0 = no jammers 

1 = one jammer 

2 = two jammers, etc. 


18-19 


Spares 


20 


Altitude in feet above which a raider will be 
considered as being at "high altitude" for pur- 
poses of obtaining a visual identification of a 
non-radiating picket. 


21 


1 -i f u p -f n i or'h -hnrJo f rnm \ ( 90 1 

raider. launching ASM at a non-radiaxing picket. 


22 


Probability of "low altitude," below X(20) in- 
put, raider launching ASM at non-radiating 
picket . 


23 


Maximum range (nautical miles) at which a raider 
is able to obtain a visual identification of a 
Blue surface unit. 


24 


Geographical plots: 

0 = no plots required 

1 = print plot of Blue positions 

2 = print plot of raider tracks 

3 = print both plots 


25-50 


Spares 


(51-60) 


(Debug switches - cleared at end of first itera- 
tion a printout will occur if input value is 1.) 


51 


Ship table. 


52 


Raid table. 


53 


SAM table. 
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X,(.. 


Inputs 


54 


Random numbers used. 


55 


Subroutine login/logout. 


56 


SAM/ ASM engagement, ASM impact, and SAM kill. 


57 


ASM detection data. 


58 


Closest point of approach and IBAT data. 


59 


SAM-ASM Intercept data (angles). 


60 


Store and Get Event- Log. 


(61-70) 


(Statistics generated by the model, ignore for 
input . ) 


61 


Total ASM's fired for all iterations. 


62 


Total ASM's failed at launch. 


63 


Total ASM's killed by SAM's. 


64 


Total ASM's failed in flight. 


65 


Total ASM's destroyed with launch platform. 


66 


Total SAM launches. 


67 


Total SAM failures. 


68 


Total SAM overkills. 


69 


Total ASM successfully launched. 


70 


Total ASM hits. 


71-74 


Spares 


75 


Number of SAM rounds remaining in Fleet. Inter 
nal use by the model. 


76 


Hazard zone radius, nm. If predicted intercept 
of target is within this range of another ship, 
then the firing ship must hold fire until the 
predicted intercept point has moved beyond this 
range . 


77 


Breakpoint altitude, ft, used to select proper 
SAM velocities, ranges, and intercept contours, 
for "high" or "low" SAM contours. 
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X£ 2 

78 

79-80 

81 

82 

83 

84 

85 



86 

87 

88 
89 

90-100 



Inp uts 

&hip position inputs coordinate system: 

0 = Polar (000 True is Y axis). 

1 = Cartesian (X = ; Y = ). 

Spares 

Radar horizon factor (default value 1.23). 

End of game time (default value 120 min). 

Launch interval for one raid, min. 

Launch interval for all raids, min. 

Number of ship for which the display of hit fre- 
quency and hit interval is desired. Enter 0 if 
analysis is not desired or if the number of 
iterations is less than 10 or greater than 100, 
in which case the histogram can not be generated 
by the model. 

Interval for initial launch of Red SSM's, min. 

Recognition delay, in minutes, for stationary 
Red raiders. For use when X(10) is not zero. 

Number of ships. 

Number of raiders (including Red aircraft which 
are only jammers, those without weapons, and Red 
surface/subsurface units). 

Spa.res 



2 . Ship Class Inputs 



The different types of BLUE ships in the simulation 
are specified by the ship class inputs. Included under one 
listing are all of the ships of a particular class (e.g. 
DDG, CG, FFG) which possess similar weapons and detection 
capabilities. A maximum of 10 different classes may be 
used. Ship class 1 is input into X(101) to X(140); ship 
class 2 into X(141) to X(180); etc. 
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X£ l 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 



116 



117 



Inputs 

- -* 

SAM/SSM battery #1 class'*" 

Sam/SSM battery # 2 class 
SAM/SSM battery #3 class 
SAM/SSM battery #4 class 

2 

Search radar # 1 class number 
Search radar #2 class number 
Search radar #3 class number 
Search radar #1 antenna height (ft) 

Search radar #2 antenna height (ft) 

Search radar #3 antenna height (ft) 

Fire control radar #1 antenna height (ft) 

Fire control radar #2 antenna height (ft) 

Fire control radar #3 antenna height (ft) 

Fire control radar #4 antenna height (ft) 

Coordination capability for this ship class: 

0 = none (override) 

1 = coordination X(14) if fleet coordination 

is in use. 

Probability of correct identification and desig- 
nation to Weapon Direction Equipment (WDE). 

Class of RTDF #2 (detection to ID and designa- 
tion to WDE). 



Enter 1 if the battery has SAM/SSM' s where class in- 
puts are listed in X(501) to X(540); enter 2 if SAM/SSM 
inputs are in X(541) to X(580); etc. 

2 

Enter 1 if the search radar class inputs for the re- 
spective search radar are contained in X(901) to X(920); 
enter 2 if the search radar class inputs are in X(921) to 
X(940); etc. Note that the Fire Control Radar class is 
listed with its associated SAM class inputs. 
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X(__i 


Inputs 


118 


Class of RTDF #3 (designation to assign FCS). 


119 


Class of RTDF #5 (firing delay). 


120 


Class of RTDF #7 (decisional delay for horizon/ 
sector scan mode -- used with X(7)). 


121 


Number of hits required to put ship out of action 


122 


Left SAM battery #1 launcher train angle limits 
(deg. true). 

♦ 


123 


Right SAM battery #1 launcher train angle limits 
(deg. true). 


124 


Left SAM battery #2 launcher train angle limits 
(deg. true). 


125 


Right SAM battery #2 launcher train angle limits 
(deg. true). 


126 


Left SAM battery #3 launcher train angle limits 
(deg. true). 


127 


Right SAM battery #3 launcher train angle limits 
(deg. true). 


128 


Left SAM battery #4 launcher train angle limits 
(deg. true). 


129 


Right SAM battery #4 launcher train angle limits 
(deg. true). 


130 


Ability to determine target at which an ASM has 
been fired: 

0 = No 

1 = Yes 


131 


Designation for non-radiating picket units using 
bistatic radar: 

0 = No; unit will be targeted using only the 

normal target list algorithm. 

1 = Yes; unit will be targeted using both 

the normal target list algorithm and the 
additional algorithm for non-radiating 
pickets . 

If some bistatic units are to be used in the 
main force, a separate ship class must be en- 
tered for them. 


132-140 


Spares 
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XL_1 


Inputs 


141-180 


Ship class 2 


181-220 


Ship class 3 


221-260 


Ship class 4 


261-300 


Ship class 5 


301-340 


Ship class 6 


341-380 


Ship class 7 


381-420 


Ship class 8 


421-460 


Ship class 9 


461-500 


Ship class 10 


3. 


SAM/SSM Class Inputs 



These inputs specify the BLUE SAM/SSM system capa- 
bilities in terms of missile flight characteristics, fire 
control radar class, salvo firing doctrine, missile relia- 
bilities and kill probabilities. A maximum of 10 different 
SAM/SSM classes may be input. SAM/SSM class 1 is input 
into X( 501 ) to X( 540) ; SAM/SSM class 2 into X(541) to X(580); 



etc . 


• 


x( ) 


Inputs 


501 


Max range, nm, against "high" altitude target. 


502 


Average speed, knots, for expected intercept 
range against "high" altitude target (converted 
by the model to nm/min). 
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X( ) 


Inputs 


530 


3 

A, nm, "high" altitude target 


504 


3 

B, nm, "high" target 


505 


N, for "high" target 3 


506 


SLOPE against "high" target 3 


507 


SEEK, nm, for "high" target 3 


508 


Max range, nm against "low" altitude target 


509 


Average speed, knots, against "low" altitude 
target (converted by the model to nm/min). 


510 


A, nm, against a "low" target 


511 


B, nm, against a "low" target 


512 


N, for "low" target 


513 


SLOPE for "low" target 


514 


SEEK, nm, for "low" target 


515 


Minimum range, nm 


516 


Maximum altitude, ft 


517 


Number of Fire Control Radar (FCR) channels for 
FCR Type 1. (0= missile not requiring a dedi- 

cated fire control system. ) 


518 


Class of FCR - Type 1 


519 


Number of FCR channels for FCR Type 2 


520 


Class of FCR - Type 2 


521 


Number of FCR channels for FCR Type 3 


522 


Class of FCR - Type 3 


3 The 
tained in 


definitions of A, B, N, SLOPE, and SEEK are con- 
the description of subroutine TIMMY in Appendix B 
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X( ) 


Inputs 


523 


Horizon sector scan capability for fire control 
radars : 

0 = No (used when X(7) set) 

1 = Yes 

input as XXX = unit digit FCR Type 1 

ten digit FCR Type 2 
hundred digit FCR Type 3 


524 


Number of SAM's per salvo 


525 


Launch cycle time, min 


526 


Probability that launcher is available (used 
only when X(2) is set to 1). 


527 


Magazine size 


528 


Spare 


529 


Probability of single-shot kill 


530 


System reliability 


531 


Reliability (launch, guidance, and fuse) 


532 


Terminal guidance time, min (if HAV. T set to 0). 


533 


Kill assessment time, min 


534-535 


Spares 


536 


Flag for Blue surface to surface missiles: 

0 = SAM 

1 = SSM 


537 


Single shot kill probability for a crossing tar 
get 


538-540 


Spares 


541-580 


SAM Class 2 


581-620 


Sam Class 3 



861-900 


SAM Class 10 
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4 . Radar Class Inputs 



Search and fire control radars for both BLUE and RED 
units are specified in this section of the inputs. Search 
radars classes are related to the ship class and raid class 
inputs. The fire control radars are used with the appro- 
priate SAM/SSM class inputs although their antenna heights 
are contained in the ship class inputs. The radars are 
specified on the basis of their detection capabilities, beam 
widths, side-lobes, radar frequency, and delay functions. 

In the case of a bistatic radar, the location and antenna 
height of the transmitter must also be specified. A maximum 
of 15 radar classes may be included. Radar class 1 is input 
into X( 901 ) to X(920) ; radar class 2 in X(921) to X(940); 



etc . 




XL..). 


Inputs 


901 


Radar Type: 

0 = Bistatic 

1 = Monostatic 


902 


x 

B, detection range in clear environment, nm/m a 
(see Appendix C). 


903 


a, detection j range with standoff jamming of 
lw/mhz, nm/ra 2 (see Appendix C). 


904 


Side- lobe ratio. Number of db that average 
side-lobe is below main beam. Enter as a nega- 
tive number. Converted by model to side-lobe 
coefficient for use in obtaining jamming effects 


905 


Maximum instrumented range, nm. 


906 


Class of reaction time distribution function: 

Detectability-to-detection for search radars 
FCS assigned to acquisition for fire control 
radars. 
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Xv 1 

907 

908 



909 

910 

911 

912 

913 

914-916 

917 



918 

919-920 

921-940 



Inputs 

Radar frequency band index. 1, 2, 3, or 4. 

Used with jamming to provide for frequency dif- 
ferences in various radars. 

Class of reaction time distribution function for 
fire control radars: 

Detectability to detection delay in horizon/ 
sector scan mode. 

Blank for search radars. 

Max Director assignment range for FCR, nm. 

Spare 

Radar beam width in degrees. (Converted to one- 
half beam width in radians by the model.) 

Spare 

Probability radar is fully operable. (Used when 
X(2 ) set. ) 

Spares 

Ship number of bistatic transmitter. Used only 
with bistatic. radars. If transmitter is not 
located on a ship, a dummy ship must be input at 
the appropriate geographical location of the 
transmitter . 

Bistatic transmitter antenna height in feet. 

Used only with bistatic radars. 

Spares 

Radar Class 2 



1181-1200 Radar Class 15 

5. Ship Position Inputs 

BLUE ship positions are input by specifying the 
class of a ship that is at a given location. The position 
may be entered as the range and bearing from force center or 
as X and Y coordinates in nautical miles. The form must be 
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the same for all ships and is dependent on the value of 
X(78). The left and right firing bearing for the ships must 
also be specified if firing sectors' are used. If a mix of 
some ships with firing sectors and others without is used, 
X(8) should be set to 1 and appropriate bearings entered for 
all ships. Ships which are not limited should be assigned a 
left bearing of 000 and a right bearing of 359.99. A maxi- 
mum of 40 ships may be used. Ship 1 is input into X(1201) 
to X(1205); ship 2 in X(1206) to X(1210); etc. 



X(, 

1201 

1202 

1203 

1204 

1205 



Inputs 

Ship class number. 

If X(78) = 0, circle in nautical miles 
= 1, X coordinate in nm 

If X(78) = 0, bearing from forces center in 
degrees true 
= 1, Y coordinate in nm 

Left sector limit, degrees true, use if X(8) = 1 

Right sector limit, degrees true, use if X(8)= 1 
In assigning sector limits to ships with no 
limitations, put 000 in X(1204) and 359.99 in 
X( 1205 ) . 



1206-1210 Ship number 2 



1396-1400 Ship number 40 

6 . Reaction Time Distribution Function Inputs 

These inputs provide appropriate delays, in minutes, 
for processing events during the simulation. A detailed 
description of the functions together with an example is 
contained in Appendix B under subroutine REACT. There are 
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seven events which may be specified by up to five different 
classes of delay functions for each event. The reaction 
functions 2, 3, 5, and 7 are related to ship class inputs. 
Functions 1, 4, and 6 are used with radar class inputs. 



# 1 : 



Detactability to Detection (5 classes) 



X( ) 


Inputs 












1401 


Delay 


in 


minutes 


for 


0.1 


cumulative 


probability. 


1402 


Delay 


in 


minutes 


for 


0.2 


cumulative 


probability . 


1403 


Delay 


in 


minutes 


for 


0.3 


cumulative 


probability . 


1404 


Delay 


in 


minutes 


for 


0.4 


cumulative 


probability. 


1405 


Delay 


in 


minutes 


for 


0.5 


cumulative 


probability. 


1406 


Delay 


in 


minutes 


for 


0.6 


cumulative 


probability . 


1407 


Delay 


in 


minutes 


for 


0.7 


cumulative 


probability . 


1408 


Delay 


in 


minutes 


for 


0 . 8 


cuiiiula Live 


probab i 1 i l y . 


1409 


Delay 


in 


minutes 


for 


0.9 


cumulative 


probability . 


1410 


Delay 


in 


minutes 


for 


1.0 


cumulative 


probability . 


1411-1420 


Class 


2 of RTDF1 











1441-1450 Class 5 of RTDF1 

#2: Detection to ID and Designation to Weapons Direc- 

tion Equipment (5 classes) 

1451-1460 Class 1 of RTDF2 



1491-1500 Class 5 of RTDF2 
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2LL_1 



Inputs 



#3: 

1501-1510 

1541-1550 

#4: 

1551-1560 

1591-1600 

#5: 

1601-1610 
A r* A -a -a r~ fs 

lO^iJL-XOJU 

# 6 : 

1651-1660 

1691-1700 

#7 

1701-1710 

1741-1750 

1751-1800 



ID and Designation to Fire Control System (FCS) 
Assignment (5 classes) 

Class 1 of RTDF3 



Class 5 of RTDF3 

FCS Assigned to Acquisition (5 classes) 
Class 1 of RTDF4 



Class 5 of RTDF4 

FCS Acquire to Fire (5 classes) 

Class 1 of RTDF5 



Class 5 of RTDF5 

Detectability to Detection for Fire Control Radar 
in Horizon/Sector Scan Mode (5 classes) 

Class 1 of RTDF6 



Class 5 of RTDF6 

Decisional Delay for Horizon/Sector Scan Mode 
(5 classes) 

Class 1 of RTDF7 



Class 5 of RTDF7 
Spares 
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7. Raid Class Inputs 



Raid classes are used to specify the launch plat- 
forms for the RED ASM/SSM and standoff jammers. The launch 
platforms can be aircraft or surface units. Submarines are 
required to be on the surface for launch and are input as a 
surface unit. The raid classes contain the offensive capa- 
bility of the RED units in terms of weapons carried, radar 
class, jamming ability, and speed. The inputs also specify 
the radar cross section of the raider and the number of hits 
to put it out of action. A maximum of 10 classes may be in- 
put. Raid class 1 is contained in X(1801) to X(1820); raid 
class 2 in X(1821) to X(1840); etc. 





Inputs 


1801 


Spare 



1802 


Number of hits to put out of action, 
for aircraft. Must be whole number. 


always 1 


1803 


Number of first type weapons carried. 




1804 


Class of first type weapon. 




1805 


Number of second type weapons carried'. 


1806 


Class of second type of weapon. 




1807 


Radar class number. 




1808 


Salvo switch (surface and subsurface 
0=1 missile per salvo 
1=2 missiles per salvo 


only) 


1809 


Launch cycle time for salvo, minutes 
subunits only. 


for surface/ 


1810 


Jamming power, w/MHz on band 1. 




1811 


Jamming power, w/MHz on band 2. 




1812 


Jamming power, w/MHz on band 3. 
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2L( 1 

1813 

1814 

1815 

1816 



Inputs 

Jamming power, w/MHz on band 4. 

Raid speed, knots (converted by the model to 
nm/min ) . 

Raid radar cross section, meters squared. 

Fourth root of raid radar cross section (computed 
by the model). 



1817 Jammer beam width in degrees. (Converted to 
one-half beam width in radians by the model.) 

1818 Jammer side-lobe ratio. Number of db that aver- 
age side-lobe is below main beam. Enter as 
negative number. (Converted by model to side- 
lobe coefficient for use in obtaining jamming 
effects. ) 



1819-1820 Spares 



1821-1840 Raid Class 2 



1981-2000 Raid Class 10 

8 . Raid Weapon Class Inputs 



The raid weapon classes specify the types of ASM/SSM 
by their flight characteristics, reliabilities, and radar 
cross section. A total of 20 different classes may be used. 
Weapon class 1 is input in X(2001) to X(2010); weapon class 
2 in X(2011) to X(2020); etc. 

X( ) Inputs 

2001 ASM/SSM speed, knots (converted by the model to 
nm/min ) . 

2002 ASM/SSM flight altitude, ft. 
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X_(_l 

2003 

2004 

2005 

2006 

2007-2008 

2009 

2010 

2011-2020 



Inputs 

Launch range"* 

Reliability (launch and- initial guidance). 

Inflight reliability (including homing and fus- 
ing). 

Radar cross-section, meters squared. 

Spares 

Weapon horizon, nm (computed by model). 

Fourth root of radar cross section (computed by 
model ) . 

Weapon Class 2 



2192-2200 Weapon Class 20 



9 . Raid Position Inputs 

These inputs contain the location, course, and alti- 
tude of an individual launch platform. If the raider has a 
jammer, the jammer target is also specified. A maximum of 
90 raiders is permitted. Raid 1 is input in X(2201) to 
X(2206 ) ; raid 2 in X(2207) to X(2212); etc. 



* (.. ) , 
2201 
2202 
2203 



Inputs 

Raid class number (1, 2, etc.) 
Range, nm (from (0,0)). 

Bearing, degrees true (from (0,0)) 



Launch Range is measured from the position (0,0). The 
user must consider this fact when dealing with dispersed 
formations and multiple target list. The distance of the 
intended launch point from force center is the proper input , 
not the range of the launch point from the primary target. 
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X( ) Inputs 

2204 Altitude, feet (less than 50 feet for surface/ 

subsurface units). 



2205 



Course, degrees true. 



2206 Jammer's jamming target: 

0 = no jam capability 

N = ship number of jamming target 
999 = force center (0,0) 



2207-2212 Raid number 2 



2734-2740 Raid number 90 

B. THE TARGET LIST ASSIGNMENTS 

The target list assignments are read in immediately fol- 
lowing the zero card indicating the end of the X vector 
inputs. Each RED raider must be assigned a target list, 
even if the raider carries no weapons. If the raider car- 
ries no weapons the target list information is not used, but' 
the input logic requires that some number be entered. There 
must be a total of three cards which are to be dimensioned 
3012. The first thirty raiders are entered with card one, 
the second thirty raiders on two, etc. For example, if the 
simulation has four RED aircraft, column 2 card 1 would con- 
tain the number of the target list for raider number 1; 
column 4 card 1 that for aircraft number 2; etc. Columns 9 
through 60 on the first card, and 1 through 60 on the next 
two cards could contain any integers desired, usually 00. 
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C. THE TARGET LISTS 

The cards that input the target lists follow immediately 
after the three cards described above. The values to be 
read in are stored in the array PLIST. This is a (40, 6) 
array in which the columns are the target lists. The array 
allows for up to 40 individual targets and six separate 
listings of targets. The target list is identified by its 
column number, and it is this number that was entered in 
NRTGT. The PLIST cards follow immediately after the three 
input cards for NRTGT mentioned above. The format for the 
PLIST cards is 10F8.4. The column entries are the total 
probability of targeting all ships up to and including the 
target ship in the given row. As an example, in the case of 

t Vvv»o<^v H*ov»rrr\-f- "1 4 c< o w H -Pi vo T3TTTTP eh i pc' f ho TDT T QT o v» o *1 rl 

j.j. x o tx x ^ w/ _i_ x. »_> L o uii x X * L- x/Xri o U o xx _i_ o t xxo x x — i kv x x u j ox kk 

appear as follows: 



1 2 3 4 5 6 



1 


. 12 


O 

o 


. 25 


0 


0 


0 


2 


.42 


. 25 


. 50 


0 


0 


0 


3 


.62 


.50 


. 50 


0 


0 


0 


4 


.82 


. 75 


1.0 


0 


0 


0 


5 


1.0 


1.0 


1.0 


0 


0 


0 


6 


0 


0 


0 


6 


0 


0 



The raiders that have been assigned target list 1 would have 
a probability of .12 of firing a missile at ship 1 each time 
it fires. It would have a probability of . 3 of firing at 
ship 2, and .20 of firing at ship 3. A raider firing a 
missile with target list 2 would have a probability of 0.0 
of firing at ship 1 and .25 of firing at any of the other 
four. A raider with list 3 would have a .25 probability of 
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firing at ship 1, a .25 probability of firing at ship 2, and 
a .50 probability of firing at ship 4. 

The input cards for the PLIST array would be structured 
as follows. Card one would contain .12 in columns 1-8, .42 

in 9-16, .62 in 17-24, .82 in -25-32, 1.0 in 33-40, and zeros 

throughout the remainder. Cards two, three, and four would 
contain zeros. Since card five starts list two, it is filled 
out in the same manner as card one, with, in this case, 
cards 6, 7, and 8 containing zeros. Since there are only 
five ships and three lists, the last card required for this 
run would be card 9. If the full 40 ships and six target 
lists were being used, a total of twenty-four cards would be 
required to input the target lists. 
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V . MODEL OUTPUT 



The model produces three types of output. The types are 
standard, special displays, and debugging. For debugging, 
there are a number of switches in the control section of the 
inputs which may be set to provide intermediary values dur- 
ing the iteration. The special displays provide one time 
outputs of geographical plots, histograms, and selected 
tables. The standard output occurs at the end of each itera- 
tion and contains an abbreviated battle history along with 
statistical data concerning the engagement. 

A . STANDARD OUTPUT 

o s*t 3.2? clstpcl output consists oi tl^o ASM t u b t o tlno SAM— 
ASM-RAID table, and statistics regarding hits, detections, 
and opportunities to engage. The contents of the tables are 
given below. The actual format of this output is provided 
in the tables of representative output at the end of this 
section . 

1. The ASM Table 

Column 1 ASM/SSM number. This number identifies 
the ASM/SSM by target number. The order 
in which the missiles were launched or 
destroyed can be determined by subtracting 
90 from the target number. 

Column 2 Status code. This code indicates the 

disposition of the missile. The code is 
as follows: 

0 missile failed on launch. 

1 missile was in flight. 
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2 missile was destroyed by SAM. The 
SAM-ASM-RAID table will indicate which 
missile ship was responsible. 

3 missile failed in flight. 

4 missile flight successful, missile 
impacted on target ship. 

5 missile destroyed when raider destroyed. 

Column 3 Time of impact. This will be the time the 
missile impacted, or would have impacted, 
unless it was lost when raider was de- 
stroyed. In this case it will be 0.0. 

Column 4 Launch raid. 



Column 5 Target ship. Target ship number unless 

the missile was lost with raider destruc- 
tion. 



Column 6 Launch time. The time of missile launch 
unless lost when raider destroyed. In 
this case it is the time of loss. 

Column 7 Weapon class. If lost with raid, 0.0. 

Column 8 X coordinate at launch. 

Column 9 Y coordinate at launch. 



2. The SAM-ASM-RAID Table 



Column 1 Target number. Targets 1 through 90 are 
raiders, 91 through 190 are missiles. 

Column 2 Launcher. The launcher code AAB, where AA 
is the ship number and B is the launcher. 

Column 3 Channel. The director-channel code AB 
where A is the director number and B is 
the channel number. 



Column 4 Time of Launch. Time of BLUE missile 
launch . 



Column 5 Kill code. The code is as follows: a 

negative number, an overkill, the number 
indicating the number of missiles in the 
salvo; a zero, a missile failure or system 
support failure; a positive number, a 
kill, the value indicating the number of 
missiles in the salvo. 
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Column 6 Time channel free. The game time that 

this launcher-channel combination will be 
available for re-engagement or reassign- 
ment . 

3. Statistical Output 

The output form of the statistics for each iteration 
are underlined. The definitions of the terms used follow 
each form. 

The ships took "H" hits : The total number of hits 

(H) sustained by the BLUE force. 

Of "X" ASMs fired, "Y" ("Z" percent) failed at launch . 
X: total number of ASM/SSM that survived to launch. 

Y: of X, the number that failed on launch. 

Z: Y/X. 

Of launched ASM, "W" ("P" percent) were killed by 
SAMs, "Q" ("R" percent) failed in flight . 

Launched ASM; from above, X-Y. 

W: the number of RED SSM/ASM shot down. 

P: W/Launched ASM. 

Q: the number of RED SSM/ASM that failed in flight. 

R: Q/Launched ASM. 

Of total ASMs available "I" ("J" percent) were 

destroyed when mother killed by SAMs . 

Total ASMs: All SSM/ASM carried by engageable raid- 

ers or raiders that reach a launch 
position . 

I: The number of RED ASM/SSM lost when raider de- 

stroyed . 

J: I/Total available. 
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”D" SAMs were launched of which "E" ("F" percent) 



failed and "G" ("H" percent) were overkills . 

D: Total BLUE SSM/SAM launched. 

E: Number of SAM/SSM evaluated as failure at inter- 

cept or system support. 

F: E/D. 

G: Overkills. Overkill occurs when the target does 

not survive to the SAM intercept point due to 
being shot down. If the target is a missile and 
the RED missile reaches its target prior to 
intercept, this intercept is scored as a failure 
and not an overkill. 

H: G/D. 

The ship statistics for each iteration are printed 
in a table labeled SHIP STATS. The format is as follows: 



Column 1 Ship number. 

Column 2 Chances to engage. The number of FED tar 
gets that will enter this ship's missile 
envelope . 

Column 3 Number of engagements. The number of 
actual missile intercept evaluations. 

Column 4 Number of targets killed. 

Column 5 Bombers detectable. Of those • targets 

that will enter the ship's missile enve- 
lope, the number that are raiders. 

Column 6 Bombers detected. Of the raiders detect- 
able, those detected. 

Column 7 ASM's detectable. Those ASM/SSM that 

will enter the ship's missile envelope. 



Column 8 ASM's detected. 

The following additional statistics are provided 
with the SHIP STATS table: 



Average detection range : The average detection 

range of each ship on those targets that 
entered its missile envelope. 
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"X" ASMs were fired on by the fleet . X is the num- 
ber of ASM/SSM that were engaged by one 
or more BLUE units. 

"X" Bombers were fired on by the fleet . X is the 
number of RED raiders that were engaged 
by one or more BLUE units. 

ASMs fired at . The number of ASM/SSM that were f irst 
engaged by each ship. 

Bombers fired at . The number of RED raiders that 
were first engaged by each ship. 

A table of the frequency of hits received by ships 
is printed at the end of the run. This table displays a 
distribution of the hit frequency over all ships for all 
runs. It constitutes a histogram of hit frequencies for all 
BLUE units. 



B. SPECIAL DISPLAY SWITCHES 



The below listed switches will provide the indicated 
output if set to one or the specified value. Setting X(2) 
to 1 will provide the equipment availability status for each 
iteration. The remaining switches provide a one-time output 
for a run. 



X(2) Indicates equipment availability check desired. 

Prints out each ship number and equipment status. 

X(24) Plot of BLUE unit positions using DISPL if set 
to 1 . 

Plot of raider tracks and launch points using 
RAIDIS if set to 2. 

Both plots are provided if set to 3. 

X(51) Prints the SHIP table at the end of the first 

iteration. The table is structured as follows: 
(see Table IV) 

. Column 1 Ship number. 
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X(52) 



Column 2 Ship class number. 



Column 3 Ship's X coordinate. 

Column 4 Ship’s Y coordinate. 



Column 5 Hits to sink. The number of RED 

SSM/ASM missile hits to put out of 
action. 



Column 6 Class address. The ship class start- 
ing index in the X array. 

Column 7 Up /down. 0 = sunk 

1 = operating normally. 

Prints the RAID table at the end of the first 
iteration in the following form: (see Table IV) 



Column 1 RAIDER The raid number, 1-90. 

Column 2 CLASS The class starting index in 

the X array. 

Column 3 TIMER (R) The time the raider was at 

the position indicated in 
columns 4 ?,nd 5 . 

The raider's X position at 
the time specified in 
column 3. 

The raider's Y position at 
the time specified in 
column 3. 

The X position of the 
raider at the time speci- 
fied as the end of the game. 



Column 4 XR 



Column 5 YR 



Column 6 XEOG 



Column 7 YEOG The Y position at the end 

of the game. 

Column 8 XL The raider X position at 

its first launch time. If 
the raider has no missiles, 
or cannot reach a launch 
position, the end of game 
position will also be indi- 
cated here. 



Column 9 YL 



The raider Y position at 
first launch. 
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Column 10 



X( 53 ) 



ALT The raider altitude in feet. 

If 50 feet or below, the 
raider is considered a sur- 
face or subsurface unit. 

Column 11 HORZ The raider radar horizon 

in miles. 



Column 12 WEP-CLASS The class number of the 

raider's longest range 
ASM/SSM. 



Column 13 WEPS 



The total number of weapons 
carried by the raider. 



Column 14 UP/DOWN The status of the raider 

at the end of the game. 

0 = shot down or sunk 

1 = operating normally. 



Prints the SAM table at the end of the first 
iteration. The table is organized as follows: 
(see Table V) 



Column 


1 


SAM 


The salvo number. 


Pr> 1 nrrm 


2 


SAM CLASS 


The SAM type class number 




- 


INDEX 


for the SAM/SSM in this 
BLUE miss-le salvo. 


Column 


3 


SHIP 

NUMBER 


The firing ship. 


Column 


4 


SAM LEFT 


The number of SAM/SSM re- 
maining in the battery's 
magazine after this salvo 
was fired. 


Column 


5 


CHI FREE 


This is the time that chan- 
nel 1 for the fire control 



radar used in this inter- 
cept will be free for re- 
engagement or reassignment 
If the missile requires 
director illumination 
throughout missile flight, 
or the missile is totally 
independent of the direc- 
tor, this is an indication 
of system availability. 
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Columns 6-10 



Channel availability times 
for systems employing a 
time sharing system of 
horning . 

Column 11 SALVO The number of missiles in 

the salvo. 

X(85) If the number of iterations is between 10 and 
100, prints a histogram of hit frequency and 
ASM/SSM arrival interval for the BLUE ship who’s 
number has been entered. 



C. THE DEBUGGING SWITCHES 

The debugging switches, X-3 , 9 , 12 , 15 , 54 , 55 , 56 , 57 , 58 , 59 , 

60, are set by assigning a value of 1 to the specified loca- 
tion in the X-array. These switches are subroutine related 
and will provide the indicated output whenever the program 
cycles through the appropriate portion of the subroutine. 
Reference should be made to the description of the subroutines 
in Appendix B when using the debugging switches. 

In the below format, the debug switch which must be set 
to 1 is listed and is followed by the subroutine name with 
its associated output. The output is in the form as it ap- 
pears in the computer printout. The parenthesis indicate 
the location of the intermediary values. 

X-3 

SUBE1 Detection range by ship ( ) on raider ( ) 
with radar ( ) in clear, jamming, and 
horizon limited environments. 

Reaction time delays for raid ( ) by ship ( ). 

Bistatic radar gain and lose contact times 
on raider ( ) by ship ( ). 
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SUBE2 



ASM/SSM failed cn launch. 



Detection range by ship ( ) on ASM/SSM ( ) 
with search radar ( ) in clear, jamming, 

and horizon limited environments. 

No detection due to reaction time delay. 

Detection range on ASM/SSM by ship ( ) on 
missile ( ) with system ( ) fire control 
radar ( ) in clear, jamming, and horizon 

limited cases. Jamming range of 999. indi- 
cates no jamming. 

Earliest search radar detection time and 
earliest fire control radar detection time 
( ) ( )• 

When missile launch time will exceed the 
launch interval set by the user, the launch 
does not take place. Output is: bomber ( ) 
could not fire as launch time was ( ) and 
all ASM's had to be launched by ( ). 



X-9 

SDBE3 EV 3 Snip ( ) target ( ) reason ( ) . 

EV 3 SAM ( ) battery ( ) channel ( ). 

PRIOR PRIOR debug data input channel ( ) output 

channel ( ) battery number ( ) input target 
( ) output target ( ). 

TIMMY Intercept data: target, firing ship, firing 
ship position, intercept position, number of 
firing zones, left and right bearings of the 
firing zones, launch bearing of SAM. 

Target, firing ship, indicator flag, reason. 

X-12 

MAIN Original X array as read into the computer. 

The target lists as read into the computer. 

The target list assignments as read into 
the computer. 
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X-15 



MAIN 

X-54 


The X-array as converted by the program. 


MAIN 


Prints out all random numbers used in the 
first iteration, in the order they were 
used. 


X-55 


Provides the computer time that the program 
enters and leaves the specified subroutine. 


EFFECT 


Login, logout 


JULIE 


Login, logout 


EVENT 


Login, logout 


GETEV 


Login, logout 


STOREV 


Login, logout 


REMOVE 


Login, logout 


TMDASM 


Login, logout 


TMDBOM 


Login, logout 


SUBE1 


Login, logout 


SUBE2 


Login, logout 


SUBE3 


Login, logout 


SUBE4 


Login, logout 


SUBE5 


Login, logout 


SUBE6 


Login, logout 


LOCTIM 


Login, logout 


PRIOR 


Login, logout 


TIMMY 


Login, logout 


RCHECK 


Login, logout 


JAMM 


Login, logout 


BIDET 


Login, logout 
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X-56 



SUBE4 


EV 4 ASM ( ) ship ( ) ASM status ( ) hits 
on ship ( ) time ( ). 


SUBE5 


EV 5 target number ( ) firing ship ( ) 
launcher ( ) salvo number ( ) time ( ). 



X-57 



SUBE1 


Detection ranges by ship ( ) on raider ( ) 
with system ( ) fire control radar ( ) in 
clear, jamming, and horizon limited cases. 

Raider ( ) surface or submarine unit, ship 
( ) has no SSM. 

Detection ranges on raider ( ) by ship ( ) 
in horizon search mode. 


SUBE2 


Reaction time delays. 

Detection ranges in horizon search mode of 
ASM/SSM ( ) by ship ( ) in clear, jamming, 
and horizon limited cases. 


TMDASM 


No engagement possible, raid ( ) target 
ship ( ) battery ( ) missile ship ( ) 
ASM/SSM ( ). 

No engagement possible, raider ( ) ship ( ) 

Engagement data - time entering and time 
leaving a given range from ship ( ). 



X-58 



JULIE 


Number of SAM/SSM in each battery. 


SUBE2 


Prints the NONO table after each raider's 
launch cycle. 

NONO table codes are contained in the de- 
scription of SUBE2 in Appendix B. 


SUBE5 


For salvo ( ) results ( ). 
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X-59 



TIMMY 


Crossing target data: target, firing ship, 
angle between target flight path and line 
of sight at launch, angle between SAM/SSM 
flight path and target flight path, angle 
between SAM/SSM flight path and line of 
sight at launch, target position at launch, 
GSINE, range to target at launch. 



X-60 



JULIE 


The events list at start of the iteration. 


EVENT 


The event retrieved. 


STOREY 


The event stored. 
The event lost. 
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NUMBER OF SAMS LEFT AT THE END OF ITERATION 1 WAS 126 
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SAM- ASM- RAID TABLE 

TGT NO. LAUNCHER CHANNEL TIME OF LAUNCH KILL CODE TIME CHANNEL FREE 

2. 21 . 17. 10.61 1. 15.36 
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10 SAMS WERE LAUNCHED OF WHICH 3 (30.00 PERCENT) FAILED AND 

0 ( 0.0 PERCENT) WERE OVERKILLS. 

TABLE tl 

REPRESENTATIVE OUTPUT FOR SAM- ASM-RAID TABLE AND THE STATISTICS 
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REPRESENTATIVE OUTPUT FOR SHIP STATS TABLE 



SHIP TABLE 
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REPRESENTATIVE OUTPUT FOR SAM TABLE 



VI. CONCLUSIONS 



The use of bistatic search radars is not a new concept. 
FADS provides a method of evaluating their potential useful- 
ness to the Navy. As demonstrated in the sample problem, 

FADS also has many applications in weapons systems analysis 
and the analysis of tactics. 

The model was structured along the lines of the FLOATS 
model in use at the Johns Hopkins University Applied Physics 
Laboratory. The additional capabilities introduced by FADS 
should provide a basis for a more realistic view of the 
present day and future tactical environments. The employ- 
ment of decoys, dispersed formations, and passive missile 
traps are but a few of the ways in which FADS might be uti- 
lized. 

Lastly, the model is ideal for use in an advanced war 
gaming course. The model is complex, but the simulation was 
written with ease of reading in mind. With rewriting the 
core requirements might be reduced, but, in its present form, 
FADS is easily read and understood, and it could be disected 
by the interested student with ease. 
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APPENDIX A 



B I STATIC RADAR SYSTEMS 



Although the bistatic radar system was described in the 
introduction as a radar system in which the transmitter and 
the receiver are not collocated it is necessary to point out 
that the separation between the transmitter and the receiver 
should.be comparable to the target distance, that is, on the 
order of several nautical miles as opposed to several feet. 
The assumed bistatic system is to provide detection and lo- 
cation information similar to that obtained with current 
search radar systems. It should in essence provide the 
user with what would be on a plan position indicator (PPI) 
scope. No provision is assumed for height finding capabil- 
ity. The basic difference between the target geometry for 
monostatic and the bistatic radar systems is shown in figure 



Al. 




Figure Ala 

Monostatic Radar Principle 



Figure Alb 

Bistatic Radar Principle 
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The bistatic radar concept is not a recent development, 
although the technology required to bring the system up to 
the capabilities described above is. The bistatic concept 
predates the monostatic radar in that the development of the 
radar duplexer was required in order to convert the original 
bistatic radar to a monostatic system. Sholnik [Ref. 10, 
p. 585] states: 

Its (the bistatic radar) principle was known 
and demonstrated many years before the development 
of practical monostatic radar. ...the first "radar" 
observations in both the United States and in 
Great Britain were made with separated CW trans- 
mitters and receivers. These early radars were 
known as wave-interference equipment but were the 
same as what would now be called bistatic radar. 

Taylor and Young of the Naval Research Laboratory 
first demonstrated bistatic radar for the detection 
of ships in 1922. Their w'ork was disclosed in a 
patent issued in 1934. The early experiments with 
wave- i nterf erence (hi static) radar led to the de- 
velopment of monostatic radar in the late 1930s in 
both this country and abroad. Further development 
was put aside after the demonstration of the more 
versatile monostatic-radar principle. Bistatic 
radar lay dormant for about fifteen years until it 
was "reinvented" in the early 1950s and received 
new interest. 

Current interest is conceivably motivated by considerations 
of electromagnetic counter-measures, defense against anti- 
radiation missiles, and deception roles. 

One method for target location with a bistatic system 
requires knowledge of the distance from the transmitter to 

the target, (R ), and to the receiver (R ), the horizontal 
x r 

angle between the base line (a line from the receiver to the 
transmitter) and the target bearing, and the elevation angle 
from the receiver to the target. However, in the simula- 
tion problem it is not necessary to determine a target's 
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location since the target's position at any time is a func- 
tion of the input parameters. It is only necessary to de- 
termine if the target is detectable and, given that it is, 
determine where along the target's track an individual radar 
is able to detect the target. 

Detection occurs when the reflected power at the re- 
ceiver is greater than the detection threshold of the re- 
ceiver for a given probability of detection. The radar 
range equation for the bistatic radar is identical to the 
well known monostatic range equation except that the power 
received varies inversely with the product of R^R 2 . Addi- 
tionally the two-way propagation losses differ since the 
path of the reflected signal is not the same as that of the 
t.ransmi tted signal. Tho range equation can be readily 
written in terms of the constants for the radar system, the 

variables R , R and the target's radar cross section, 
x r 

For comparison purposes both the monostatic and the bi- 
static detection range equations are given below: 



Monostatic P 



P t G 2 A 2 a 



r (4 tt) 3 L 2 L R 4 
v ' p s 



Bistatic P 




where : 



R = range to target (monostatic) 



R 



= range from transmitter to target (bistatic) 



x 



R 



= range from receiver to target (bistatic) 



r 
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= power of transmitter 
P^, ' = power received at the receiver 

G = gain of monostatic antenna 

= gain of bistatic transmitter antenna 
G^ = gain of bistatic receiver antenna 

X = radar wavelength 

o = target's radar cross section (see below for dis- 

cussion of bistatic and monostatic approximation) 

L = one way propagation losses for monostatic 

P 

L (t) = propagation losses from transmitter to target 
^ (bistatic) 

L (r) = propagation losses from target to receiver (bi- 



P 



static ) 

= system losses 



By defining as the minimum power required by the re- 
ceiver in order to- obtain a given probability of detection 
with typical propagation and system losses and by grouping 
constants of the radar system, the above equations can be 

i 

simplified to read: R = (3a 4 and R^R 2 = 3 4 cf where $ is the 

range at which a particular radar system can detect a one 
square meter target in a non- jamming environment. Values of 
3 for current U.S. Navy radars are listed in Ref. 9. 

When standoff electronic noise jamming is present, the 
range equations are defined in a similar manner: 



Monostatic : 



Bistatic 



P.G a 

n 4 _ t 

bt 4ttB (S/J)EJP 



x x 

P.G 2 G , 2 o 

R 2 R 2 = • - — 

bt x 4 ttB (S/J)ZJP 
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where : 



R b '= burn through range rnonostatic 

B = receiver band width 
r 

S/J = signal-to- jamming ratio (includes propagation 
losses) 

EJP = sum of jamming power arriving at the receiver 

The jamming level from an individual jammer is obtained 
by multiplying the power of the jammer by the appropriate 
coefficient for the geometry of the radar and jammer main to 
side-lobe ratios and dividing by the square of the distances 
between the jammer and the radar. The individual jamming 
levels are then added to obtain the sum of the jamming power 
arriving at the radar. 

As in the case with the non- jamming environment, the 
above equations can be simplified by defining a term a which 
contains the constants of the radar system and which is 
equal to the range at which a one square meter target can be 
detected in the presence of one watt of jamming power per 
megahertz. The monostatic range equation in jamming is: 

R bt = a^(cr/ZJP)* 
and the bistatic equation is 

R x R r = a( 0 /£ Jp )^ • 

In solving the bistatic detection problem in the simula- 
tion algorithm, R and R have only the target's position in 

x r 

common. It is therefore necessary to state the product R ,R 

X 27 

as a function of the target's position and solve the 
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detection equation in terms of the time at which the target 
can be detected. If a constant CHAR is defined to be the 
product of 3 1 * times a for a clear environment and a 2 times 
a divided by the sum of the jamming power for the jamming 
situation, the detection equation then becomes R^R^; = CHAR. 
The radar is able to detect a target whenever the range 
product squared is less than or equal to the constant CHAR. 
By expressing R and R in terms of time, the solution of 
the bistatic equation R 2 (t)R 2 (t) -- CHAR = 0 provides the 

X T 

times at which detection can occur and then be lost as the 
target moves along its track relative to the transmitter and 
the receiver. The below figure is a layout of the track and 
detection points that is illustrative of the detection equa- 
te on : 




Figure A2 

Bistatic Detection Geometry 
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Labels are as follows: 



Receiver 


Transmitter 


Value 




XR 


XX 


X Coord 


of position 


YR 


YX 


Y Coord 


of position 


ZRN 


ZXN 


antenna 


height in nautical mill 


XRC 


xxc 


X Coord 


of track CPA 


YRC 


YXC 


Y Coord 


of track CPA 


CPR 


CPX 


Distance 


to CPA 


Track 








Xi,Yl 




initial 


position 


Z1 




altitude 


in nautical miles 


XO, YO 




terminal 


position 


V 




velocity 


in nm/min 


Vx 




veloc.j tv 


in X direction 


Vy 




velocity 


in Y direction 


XD 




X Coord of detection position 
= (Xl+VxT) where T is the time 
of detection 


YD 




Y Coord of detection position 
= (Yl+VyT) 



The equation R 2 (t)R 2 (t) - CHAR = 0 is expressed in terms 

X 27 

of the stated positions; time is the independent variable. 
For example: 

R 2 (t) = (XR-XRC ) 2 + ( YR-YRC ) 2 + (ZRN-Z1 ) 2 + (XRC- (Xl-i-VxT ) ) 2 
+( YRC- (Yl+VyT ) ) 2 

= (Vx 2 +Vy 2 )T 2 -2( Vx(XRC-Xl )+Vy ( YRC-Y1 ) )T+ (CPR) 2 
+ ( XRC-X1 ) 2 + ( YRC-Y1 ) 2 

= at 2 +dt+e 
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In a similar manner R 2 (t) = AT 2 +BT+C where B = -2(Vx(XXC-Xl ) 

X . 

+Vy ( YXC-Y1 ) )T and C = (CPX ) 2 + (XXC-X1 ) 2 + ( YXC-Y1 ) 2 . 

The range equation in terms of t then reduces to 



T 4 + (B+D) T 3 + (C+ e 



A 



+ 



BD )t , +( BE 



-CD 



A 2 



)T+ 



CE - CHAR 
A 2 



0 



The real roots of this fourth order polynomial in ascend- 
ing order are the times at which the target is first detect- 
able, lost, redetected and then lost again. In the case 
where there are no real roots, the target is not detectable 
along its track. If only one pair of real roots exists, the 
radar has only one detection zone against the target. 

The multiple detection zones occur when the target pene- 

trates the contour of the equation R R = CHAR“ in more than 

1 x r 

two locations. These contours plot as ovals of Cassini in 
the plane containing the base line between the transmitter 
and the receiver. The contours, shown in Figure A3, repre- 
sent curves of constant sensitivity for the bistatic radar 
system and a given target. In the case of low sensitivity, 
the curves consist of two disjunct ellipses with the trans- 
mitter and receiver at the furthest foci. As sensitivity is 
increased the curves become ovals about the receiver and the 
transmitter and eventually become an ellipse with the trans- 
mitter and receiver at the foci. If this ellipse is rotated 
about its major axis, a prolate spheroid is generated which 
defines the detection surface for the bistatic radar. 
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Figure A3 

Curves of Constant Sensitivity for Bistatic Radar System 

In addition to the requirement that a solution to the 
radar range equation exists, the target must, also be above 
both the transmitter and the receiver radar horizon for a 
detection to occur. The radar horizon restrictions can be 
seen in figure A4, which uses the same track as presented in 
figure A2, along with an additional axis representing the 
aircraft position with respect to time. 

The positions in figure A4 are defined in figure A2. 

HORX and HORR are the distances to the transmitter and re- 
ceiver radar horizons for their respective antenna heights 
and the target's altitude. TR1 is the time the target first 
enters the radar horizon of the receiver; TR2 is the time it 
departs. TX1 and TX2 are similar times for the transmitter. 
The time that the target is simultaneously above both hori- 
zons is the interval (T®,T®) where T® = max (TR1,TX1) and 
T® = min (TR2,TX2). 
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The intersection of the joint horizon time interval 
(T°,T°) with the time intervals obtained from the range de- 
tection equation provides the time of detection and the time 
of lost detection for a target. 

These times, in conjunction with the ranges which are 
related to the times, are then used in the simulation to de- 
termine "shoot events" for the defensive forces and to accu- 
mulate statistics on the detection ranges. 
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APPENDIX B 



This 


DESCRIPTION OF THE PROGRAM 
appendix provides a description of the FADS com- 



puter program. The MAIN portion of the program is presented 
first. The principal subroutines follow with the minor sub- 
routines and functions covered last. 

1. OVERVIEW 

This paragraph gives an overview of the functions of 
the individual sections of the program. 



MAIN 


The MAIN program reads in the input values, com- 
piles statistics and controls the simulation through 
the specified number of iterations. 


EFFECT - 


Subroutine EFFECT controls the individual itera- 
tion by selecting the proper subroutine to play 
the next' event in the events list. 


JULIE 


Subroutine JULIE initializes or resets values for 
each iteration. 


EVENT 


Subroutine EVENT controls the input and retrieval 
of events from the events list. 


REMOVE - 


Subroutine REMOVE retrieves higher priority tar- 
gets from the events list. 


TARGET - 


Subroutine TARGET returns all of the "shoot" events 
which occur at a specified time. 


SUBE1 


Subroutine SUBE1 computes detections of the raiders 
by the BLUE forces and determines "shoot" events 
against the raiders. 


SUBE2 


Subroutine SUBE2 computes the RED ASM/SSM launches 
and subsequent detections followed by the "shoot" 
events against the ASM/SSMs. 


SUBE3 


Subroutine SUBE3 evaluates the engageability of a 
target based on the situation at the proposed BLUE 
missile launch time. 
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SUBE4 

SUBE5 

SUBE6 

TIMMY 

JAMM 

BIDET 

DTIMD 

CPA 

REACT 

PRIOR 

RCHECIC 



TMDASM 



TMDBOM 



ZLNKA 

TJSORT 



Subroutine SUBE4 determines ASM/SSM impacts on 
BLUE ships and evaluates out of action status for 
them. 

Subroutine SUBE5 evaluates the BLUE missile inter- 
cepts. 

Subroutine SUBE6 releases the fire control radar 
at the time of intercept. 

Subroutine TIMMY computes the predicted intercept 
point and flight time for the BLUE SAM/SSM and 
determines whether the intercept is within the 
missile's capability. 

Subroutine JAMM determines the range at which a 
monostatic radar can detect a target in the pres- 
ence of jamming. 

Subroutine BIDET computes detection ranges and 
times of detection for bistatic radars in a clear 
or jamming environment. 

Subroutine DTIMD determines the real and imaginary 
roots of the fourth order polynomial used to solve 
for the bistatic radar detection times. 

Subroutine CPA computes the coordinates of the 
closest point of approach of a target to a unit. 

Subroutine REACT computes reaction time delays. 

Subroutine PRIOR determines the most threatening 
target among those that a defensive system may en- 
gage at a given time. 

Subroutine RCHECK evaluates a prospective target 
ship to determine if the target is above the 
raider's radar horizon and within the raider's 
maximum theoretical radar range. 

Subroutine TMDASM determines if and when a RED 
target will close within a specified range of a 
BLUE unit. 

Subroutine TMDBOM determines the time when a RED 
raider will come within a given range of force 
center . 

Subroutine ZLNKA is a sorting routine. 

Subroutine TJSORT is a sorting routine. 
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LOCTIM - Subroutine LOCTIM computes the raider position at 
a given time. 

ASMTIM - Subroutine ASMTIM computes the ASM/SSM position at 
a given time. 

RANDOM - Subroutine RANDOM generates uniform (0,1) random 
numbers . 

IIISTG - Subroutine IIISTG provides a histogram and asso- 
ciated statistics of hits, hit intervals, and hits 
per minute for a specified BLUE unit. 

DISPL - Subroutine DISPL provides a geographical plot of 
the BLUE unit positions. 

RADIAN - Function RADIAN converts angles in degrees to 
angles in radians. 

RAIDIS - Subroutine RAIDIS provides a geographical plot of 
the RED raider tracks and launch points. 

SET I MX - Subroutine SETIMX provides the internal computer 
clock time. 

2. MAIN 

The simulation is initialized by MAIN which also com- 
piles iteration results for the summaries. The MAIN program 
begins the first of its three principle functions by zeroing 
specified arrays. It then proceeds to read in the X array, 
the raider target list assignments (NRTGT), and the target 
lists (PLIST). The second principle function is to enter 
those initial values that will remain constant through all 
iterations into the various tables. MAIN's third principle 
function is the accumulation and printing of the statistical 
results . 

The initial step in performing the above functions is 
the zeroing of the X array. In the next step, data are read 
into the X array by reading data cards in the format 
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( 110 , 5F10 . 3 ) . The information is first read into an integer, 
IRT, and the array RXATA(5). From there, using IRT as the 
initial X address, MAIN fills the next five X array slots 
sequentially. X array addresses that are to remain 0.0 need 
not be read in unless they are within the neighborhood of 
the five addresses contained on an input card. In this case 
they must be specified 0.0. The sequence of reading data 
cards and filling the X array will continue until a card is 
read with the value 0 for IRT and 0.0 for the remainder of 
the card input. This is the flag for MAIN to proceed to the 
next step. 

The program then obtains from the values previously 
read into the X array the starting values for ISEED, the 
seed value for the random number generator, and NI , the num- 
ber of iterations to be played. End of game time is deter- 
mined and the X array is printed as it wa.s read in. MAIN 
now computes the weapon (RED ASM/SSM) horizon as follows: 

HORZ = (CURVATURE FACTOR) x (CRUISING ALT)^ 

HORZ is in miles and ALT is in feet. MAIN next converts all 
courses and bearings from degrees true to degrees in the 
cartesian coordinate system using the relationship: 

DEG(C) = 450. -DEG(T) 

All velocities are converted from knots to nautical miles 
per minute. Jammer and radar beam widths are converted to 
half beam widths in radians, and MAIN finds the fourth root 
of the raider and RED weapon radar cross sections. Next the 
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target ] 1st assignments and target lists are read in, and, 
if X(12) is set, MAIN prints the X-array, PLIST, and NTRGT as 
read in. 

MAIN next fills out the ship table and the raid table. 
It counts up the number of BLUE missiles (both SAM and SSM) 
available, and, if X(15) has been set, MAIN prints the con- 
verted X array. If X(55) is set, MAIN calls SETIMX to start 
the internal timer. 

MAIN now ca.lls EFFECT to control the first iteration. 
Upon completion of the first iteration, MAIN prints the ASM 
table and the SAM-ASM- RAIDER table. It then computes and 
prints various statistics for this iteration and computes 
and prints cumulative statistics for all iterations thus far 



comp 1 si ec 



T -P 



V/C1 \ 
v. '-'- L J 



gHTP tnhi 



the end of iteration one. If X(52) is set the RAID table is 
printed and if X(53) is set the SAM table is printed. If 
X(54) is set, all the random numbers used in the first iter- 
ation are printed, in the order of use. 

At the end of iteration one, MAIN sets flags 9, 52, 53, 
55, 56, 57, 58, 59, and 60 to 0. MAIN then loops back and 
calls EFFECT to begin the next iteration. Upon completion 
of the last iteration, additional statistics are printed 
and the program stops. 



3. SUBROUTINE EFFECT 

EFFECT controls each iteration. If X(55) is set for 
debugging, EFFECT calls subroutine GETIMX on entry and exit 
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to login and logout. EFFECT next calls JULIE to set the 
tables and arrays that were altered during the previous 
iteration back to their proper starting values. The simula- 
tion begins when EFFECT calls subroutine GETEV for the first 
event to be played. GETEV supplies EFFECT with the type of 
event to be played, NEXTEV, the time of the event, and a 
code identifying the type of event. If the event type re- 
turned by GETEV is 0, the iteration is over and the sub- 
routine exits. If the event type is greater than 6, the 
program prints an error message and terminates. Otherwise 
EFFECT calls one of the following subroutines to continue 
the play: 

SUBE1 for a RED raider detection and intercept 
event . 

SUBE2 for a RED missile launch-target-detect- 
intercept event. 

SUBE3 for a BLUE missile firing event. 

SUBE4 for a RED ASM/SSM impact event. 

SUBE5 for a BLUE missile kill event. 

SUBE6 for a "set fire control channel free" 

event . 

Upon return from one of the SUBE's, EFFECT again calls GETEV 
for the next event to be played. Event types used elsewhere 
in this manuscript correspond to the SUBE subroutine which 
plays the event. 

4. SUBROUTINE JULIE 

Subroutine JULIE initializes or resets those items that 
are required for each iteration. If X(55) is set, JULIE 
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calls GETIM upon entry and exit. JULIE is called by EFFECT 
and calls RANDOM, TMDBOM , STOREV, and EVENT. 

EVENT is called to zero the events list, and to enter 
the end of game event into the events list. The following 
counters are then set to zero; NASM, NSAM, IFAG, IFBG, ITEMS, 
IFLASK, and IFIR3T. The following arrays are zeroed; JAM, 
IHIT, ISAG, ISBG, STATS, TGTAV, IGAGE, MONO, JACK, and 
SAMASM. The RAID table and SHIP table are reset and the SAM 
table is zeroed. The missile batteries are reloaded and, if 
X(2) is set, the availability of each BLUE radar a.nd launcher 
is determined. JULIE does this as follows. The first check 
is on the search radars; if at least one is operating, JULIE 
determines which fire control radars are operable. If no 
search rada.r is available, the program checks to see whether 
the horizon search mode for fire control radars is author- 
ized, X(7), and if so, whether the ship in question has 
radars with this capability. If the ship has a search radar 
and fire control radar combination, or a search radar com- 
bined with active, HOJ, ARM, or some other type missile 
which does not require a fire control radar, the launch 
availability is checked. If the search radars are inopera- 
tive, the fire control radars are horizon search capable, 
and at least one is operable, the launchers are checked for 
availability. If the ship has no means of firing or of de- 
tection, it is marked a non-combatant and the next unit is 
evaluated. If the launcher can be coupled with a search/ 
detection/acquisition system, the ship is marked as a G-ship, 
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the equipment status is entered in the NADAR array, and the 
next unit is evaluated. 

For each RED raider, JULIE calls STOREV to store a 
type 1 event at time 0. , and, if the raider carries ASMs or 
SSMs, a type 2 event is stored for launch time as determined 
by TMDBOM . 

Finally, JULIE totals the number of BLUE SAM(SSM) bat- 
teries available, sets the RELOAD array to zero, and, if 
X(60) is set, prints the events list. 

5. SUBROUTINE EVENT 

Subroutine EVENT is called by JULIE. It enters the 
end of game time in the events list. If X(55) is set, 

EVENT calls GETIMX upon entry and exit. The event counter 
LOVER is set to one and the events list is zeroed. End of 
game time is set equal to the integer IT using the equation: 

IT = 1000 x (T + 0.0005) 

where T is equal to the value stored in X(82). The end of 
game event, a type 0 event, is then entered in the events 
list along with the time IT. EVENT then returns to JULIE. 

a. Entry GETEV 

GETEV is an entry into Subroutine EVENT. It is 
called by EFFECT and calls GETIMX at entry and exit if X(55) 
is set. GETEV removes the first entry in the events array, 
IVENTS (1,1) and IVENTS (1,2), converts these back into 
event time, type, and participating units, and sends these 
three pieces of information back to EFFECT. The counter 
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LOVER is decremented and the entire events list is moved up 
one position. If the debugging switch X(GO) is set, GETEV 
prints out TIME, NEXTEV, and IUNIT. 
b. Entry STOREY 

STOREY is an entry in Subroutine EVENT and its pur- 
pose is to place the given event in its proper place in the 
events list. STOREV is called by JULIE, SUBE1, SUBE2 , SUBE3 , 
SUBE5, and PRIOR. If X(55) is set, STOREV calls GET I MX upon 
entry and exit and prints the login/logout time. The inputs 
to STOREV are the time of the event, the type of event, and 
the code for the units involved. STOREV first checks to see 
whether the time the event is to be played exceeds the 
allotted game time. If it does, STOREV prints a message 
giving the time of the last event, type of event, and player 
code. If LOVER is greater than or equal to 1000, the event 
is lost and STOREV prints a similar message. If both tests 
are passed, STOREV combines the event type and player code 
into one 9 digit integer (ICODE). Time is converted to an 
integer using the same equation as is in EVENT. STOREV 
places the event in the events array by time. If the time 
is not unique, it proceeds to compare the stored values of 
ICODE to the one to be entered. These are ordered from 
smallest to largest. When the proper place is found, LOVER 
is incremented, and the remainder of the array is moved up 
one address. If no place is found to insert the event, in- 
cluding even just right before the end of game event, an 
error message is printed and the game stops. If the debug 
switch X(60) is set, STOREV prints the event stored. 
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6. SUBROUTINE REMOVE 

Subroutine REMOVE is called by PRIOR when PRIOR dis- 
covers a more threatening target for the missile system 
being evaluated. If X(55) is set, REMOVE calls GETIMX upon 
entry and exit. REMOVE receives the code IUNYT from PRIOR. 
IUNYT is the (event type/player) combination that is stored 
in the events array. IUNYT was found by TARGET and evalu- 
ated by PRIOR as the next event to be played. REMOVE 
searches the events array for a matching code group. If it 
finds the group it removes both the code and the time and 
then moves the remainder of the array up one, and exits. If 
REMOVE can not find the code, it prints an error message. 
REMOVE then prints the events array and exits. 

7. SUBROUTINE TARGET 

Subroutine TARGET is called by PRIOR. Its purpose is 
to provide PRIOR with all Type 3 events which occur at the 
present game time. The input for TARGET is time and the 
outputs are IUNIT and MORE. If IFIRST is equal to zero, 
TARGET sets it equal to one. TARGET then converts time to 
the integer format of EVENT and compares it to the time in 
the array location indicated by IFIRST. If the times are 
equal, TARGET checks the event code. If it is a Type 3 
event, TARGET increments IFIRST, sets MORE equal to one, and 
exits with the event code. If the times are not equal, or 
the code is not Type 3, MORE and IFIRST are set to zero and 
the subroutine exits. If MORE was set to one, PRIOR will 
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call again until MORE is set to zero. If X(55) is set, 

TARGET calls GETIMX each time it is entered or is exited. 

8. SUBROUTINE SUBE1 

SUBE1 is the RAIDER detection subroutine. It is called 
by EFFECT and in turn calls subroutines CPA, JAMM, STOREV, 
REACT, LOCTIM, RANDOM, TMDASM , BIDET, and, if X(55) is set, 
GETIMX on entry and exit. SUBE1 tracks the raider from its 
start of the game position to its end of game position and 
determines which BLUE units, if any, will be able to engage 
the raider. If it is determined an engagement is possible, 
detection delays, command and control delays, and equipment 
availability are checked to determine whether or not to 
enter a Type 3, or shoot, event in the events list array. 

SUBE1 first sets flags NASA, NSURF, and N J . If X(16) 
is set, SUBE1 checks the ASM counter to see if the BLUE 
force has been fired upon. If X(16) is set and NASM is zero, 
SUBE1 stores a Type 1 event at game time plus 2 minutes and 
then exits. Otherwise, the raider's altitude is checked to 
determine whether the raider is an aircraft or surface/sub- 
surface unit, and then the applicable raider characteristics 
are drawn from the X- array. 

Each BLUE ship is next queried to determine whether it 
can engage the raider. If the ship has been sunk, or if it 
is a non-missile ship, the NONO table is set to 7 and the 
next unit is evaluated. The maximum missile range of the 
operable systems on the BLUE unit is obtained. The dis- 
tinction between surface and air targets is made by comparing 
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target altitude to the fifty-foot ceiling for enemy surface 
units and the fifty-foot floor for enemy aircraft. A mis- 
sile system must be designated as a surface to surface, or 
surface to air system when the system characteristics are 
input. There is no cross matching of systems to targets. • 
SAM's are fired only at aircraft and ASM's, SSM's only at 
ships . 

TMDASM is called to determine whether the target will 
enter the BLUE unit's missile envelope. If not, the NONO 
table is set- to 9 and the next unit evaluated. If the tar- 
get will enter the envelope after the end of game time, the 
NONO table is set to 8. 

When it has been determined that the target can be en- 
gaged, the target CPA is obtained and the ship's engagement 
opportunity counter is incremented. 

The detection phase begins with a check of the search 
radars. Each BLUE unit is allowed up to three. If the 
radar is opera.ble, the clear detection range, jamming envi- 
ronment detection range, combined radar horizon, and theo- 
retical maximum radar detection range are compared. The 
minimum of these is the radar's detection range for the tar- 
get in question. This range is compared to that of the other 
search radars on board, and the maximum is taken as the 
ship's detection range. 

The between radar comparison is made by converting de- 
tection ranges into time and adding in the reaction time 
delays obtained from REACT. 
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If there is no search radar operational, SUBE1 checks 
to see whether the target can be detected in the horizon 
search mode. If there has been no detection and the BLUE 
unit is operating a bistatic search radar, the subroutine 
checks to see whether BIDET indicated two detection zones. 

If so, the search radar portion of SUBE1 is redone using the 
second bistatic zone. If there remains no detection, or no 
radar capable of detection, the NONO table is set to 2 and 
the next BLUE unit evaluated. 

Given a search radar detection, SUBE1 proceeds to eval- 
uate the missile systems for fire control acquisition. This 
is done by obtaining the maximum fire control radar range 
from among those radars available, and converting this to a 
time as was done in the search radar algorithm. If the mis- 
sile of a system requires no fire control radar, the maximum 
missile range is used as that system's acquisition range. 

Having determined the maximum acquisition range, the 
probability of correct identification is compared to a ran- 
dom number drawn from RANDOM. If the ID check declares the 
target non-hostile, the NONO table is set equal to 3 and the 
next BLUE unit is evaluated. Otherwise, the time the target 
will be acquired is determined by calling TMDASM. 

The reaction times for weapons direction equipment de- 
lays, target assignment, and acquisition (if the missile re- 
quires a fire control radar) are obtained from REACT and 
added to the search radar detection time. The last time to 
shoot is determined using the time the target leaves the 
maximum fire control range, and the SAM/SSM velocity. 
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The firing delay is obtained from REACT and added to 
the greater of the detection time and the acquisition time. 
The "shoot" time is compared to the" last time to fire, and 
if the shot is still possible, the program stores a shoot, 
Type 3, event for each director- launcher combination in op- 
eration. If the missile requires no director, the number 7 
is stored as the channel number, and only one shoot event is 
stored for that system. If the unit requires a director, no 
channel number is assigned. 

If the delays have precluded a firing, SUBE1 checks to 
see whether there was a second bistatic detection. If so, 
it returns to the search radar portion of the algorithm and 
proceeds as above evaluating the second bistatic detection 
zone. Otherwise it sets the NONO table equal to 5 and proc- 
esses the next BLUE unit. 

The horizon search portion of SUBE1 is entered as de- 
scribed above. The program checks the NADAR table to see if 
the directors can operate in the horizon search mode. If 
so, the program proceeds in the same way as the search al- 
gorithm, using the fire control radar parameters. Having 
determined the detection time, the decisional and reaction 
delays for horizon scan mode are obtained from REACT. The 
shoot time is determined to be the detection time plus the 
acquisition time, plus the maximum of the firing time and 
decision to fire time. A range check is made by calling 
TMDASM and if the time to fire is less than the last time to 
fire, a Type 3 event is stored for each director operable in 
the horizon search mode. 
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After all .the BLUE units have been processed, SUBE1 
exits calling GETIMX if flag X(55) has been set. 

9. SUBROUTINE SUBE2 

SUBE2 is the ASM/SSM launch and detection subroutine. 

It is called by EFFECT and in turn calls RCHECK, REACT, 

STOREV, BIDET, RANDOM, TMDASM, TJSORT, CPA, LOCTIM, JAMM, 
and if X(55) is set, GETIMX on entry and exit. 

The inputs to SUBE2 are launch time and raider number. . 
The subroutine first checks the RAID table to see whether 
the raider has already been shot down or sunk. If so, the 
subroutine exits. Next the program checks the flag IFLASK 
to determine whether any RED ASM/SSM' s have been fired. If 
IFLASK is greater than zero, the subroutine proceeds to 
check whether the launch time is greater than TMAX and exits 
if it is. If IFLASK equals zero, SUBE2 sets IFLASK to one 
and TMAX, the time by which all ASM's must be launched, to 
present time plus X(84). SUBE2 now retrieves the raider 
speed and, if the raider is an aircraft, determines the ASM 
launch interval, sets IFLASK equal to 2, and determines the 
total number of ASM's on board. If the raider is a surface 
unit; and no ASM's have been fired, IFLASK is set to 0. 

SUBE2 determines which of the two possible missile types 
has the longest firing range, (firing range is measured from 
coordinate (0,0)). The missile parameters are drawn from the 
X vector, and the launching sequence begins. 

The launching of ASM's is done by drawing a series of 
uniform (0,1) random numbers, depending upon how many missiles 
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are to be launched, and ordering them. The exact launch 
time of each ASM becomes the present time, plus the next 
random number times the launching interval. The launching 
of SSM's from surface units begins with the first launch at 
time plus some portion of the interval specified in X(86). ' 
Subsequent surface launches are at fixed increments, X(1809). 

Having determined when, and in the case of surface to 
surface launches, how many RED missiles are to be launched, 
the targeting phase beings. LOCTIM is called to insert the 
raider position at launch into the RAID table. The raider 
target list is located and a temporary target list is con- 
structed. Each ship on the raider's target list is cycled 
through RCHECK to determine whether it is afloat, above the 
radar horizon and within radar range. If all these tests 
are passed, the ship is entered on the temporary target list 
by putting the values from the original target list in the 
next slot of the temporary target list. Once the target 
list has been formed, a random number, Z, is drawn from 
RANDOM and the first ship with a cumulative probability of 
being targeted greater than Z is designated the target. If 
there is none, the first ship on the list is taken. A check 
is made to see whether this ship has been designated a non- 
radiating picket. If not, SUBE2 proceeds to the detection 
phase. If the target chosen is a non-radiating picket, 
checks are made to see whether the raider can correctly 
identify the target. This is doen by determining whether 
the raider is within X(23) miles of the target. This 
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measurement is done in the horizontal plane. If the raider 
is close 'enough to the target, RANDOM is called for a random 
number to compare with the probability of correct identifi- 
cation. If the picket fails to qualify as a target, it is 
eliminated from the target list and the process starts again. 
If, at any time the situation arises where there are no tar- 
gets on a raider's target list, launching stops and an event 
Type 2 is stored for time plus three minutes. 

Upon acquiring a target, the ASM/SSM is considered 
launched. A random number is drawn and compared to the 
launch reliability. The ASM table is filled out and, if 
launch was unsuccessful, the NONO table is set equal to 1 
and the next missile processed. Upon an unsuccessful launch, 
a check is made to see whether there is to be a second SSM 
in the same salvo. If so, the time is not advanced for 
this launch. 

Assuming a successful launch, the detection phase be- 
gins. The target ship's position is stored in the RAID 
table for use in BIDET, but otherwise, the detection se- 
quence follows that outlined in SUBE1. 

When all of the first type of missiles have been launched, 
and if the raider is carrying a second type of missile, 

SUBE2 decides whether or not the raider has reached its sec- 
ond launch point. If so, the appropriate parameters are ex- 
tracted and the launching continues; otherwise an event Type 
2 is stored for the new launch time and the subroutine exits. 
When there are no second type missiles, the subroutine exits 
after the last launch. 
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The exit, sequence consists of checking X(5S), (if set 
the NONO' table is printed), and X(55), if set, (the exit 
clock time is printed). The NONO table codes for SUBE1 and 
SUBE2 are: 

0 System can engage the target. 

1 The ASM failed at launch. 

2 The ship did not detect this target. 

3 The missile or raider was not identified as a 
target suitable for engagement. 

4 The ship can not acquire the target. 

5 Reaction delay was too long. 

6 The missile will impact prior to coming into 
max SAM range. 

7 Non missile ship. 

8 Detection will not occur until after end of game. 

9 Can not engage. 

10. SUBROUTINE SUBE3 

This subroutine evaluates the engageability of the tar- 
get in light of the current situation at the proposed launch 
time. SUBE3 is called by EFFECT, and it in turn calls 
PRIOR, STOREV, CPA, TIMMY, and, if X(55) is set, GET I MX on 
entry and exit. 

The inputs to SUBE3 are event time and a participant 
code. The code is an integer of the following type, IIIJJKLM, 
where III is the target number, JJ is the firing ship, K is 
the missile system, L the director, and M is the channel- 
number. The subroutine breaks down the code and then checks 
the appropriate portion of the SHIP table to determine 
whether the firing ship has been sunk, (if X(5) has been 
set, this check will indicate the ship lias not been sunk 



110 



even if it has absorbed more hits tha.n are considered nec- 
essary to put it out of action). If the ship has expended 
the magazine of this missile system, the subroutine exits. 
Having established the credibility of the firing ship, SUBE3 
determines whether the raider has been shot down or sunk, • 
and if so, it exits. 

SUBE3 now determines the maximum missile range for the 
system under consideration and compares this with the raider's 
CPA. If the target will not enter this maximum missile en- 
velope the subroutine exits. PRIOR is called to determine 
whether this ship/system/director combination has another 
shoot event at the present game time and, if so, if that tar- 
get is more threatening. PRIOR logic will consider the 
closest target as most threatening. If PRIOR returns the 
same target, the subroutine continues. If not, it checks 
out the target as above and exits or continues as appro- 
priate. 

If the target has not been destroyed, SUBE3 checks the 
coordination set-up. If both the fleet and the firing ship 
are coordination capable, the subroutine checks the target 
availability. If the target is under engagement, the sub- 
routine calls STOREV to store a Type 3 event at the time the 
target will become available. If the target remains engage- 
able, SUBE3 checks the NONO table. If the value 0 is found 
the event continues. If not, the subroutine exits. 

The subroutine now goes to the channel assignment phase. 
If the SAW or BLUE SSM requires no director, channel number 
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7 has already been assigned in SUBEi or SUBE2, and the flag 
ICOT is sfet to 0. If the missile requires a director for 
part or all of its flight, and a channel has been assigned 
on a previous run through SUBE3, the subroutine checks to 
see whether the assigned channel is busy. If it is, ICOT is 
set to 1; if not, ICOT is set to 0. If the missile has not 
been assigned a channel previously, the availability of a 
channel is checked. A free channel, if found, is assigned 
and ICOT is set to 0; otherwise a search is made to deter- 
mine which channel will first become available and this one 
is assigned. If at this point all channels are busy and the 
director is required for the full SAM/SSM flight, STOREV is 
called and a Type 3 event stored for the time the director 
will become available. 

If the systems all indicate they are engaged in an 
intercept and upon searching the SAMASM and SAM tables no 
record of this can be found, the program terminates after 
printing an error message. 

SUBE3 now proceeds to check whether the firing ship 
itself is engaging the target. If it is the subroutine exits. 
If the target is an ASM or SSM, the subroutine checks to see 
that the target will not impact prior to BLUE missile launch. 
If it will, the program exits. 

At this point, if ICOT equals 1, the channels were all 
considered busy and the BLUE missile only requires the di- 
rector for terminal homing. SUBE3 checks the SAMASM listing 
to determine when the channel will be available and stores a 
Type 3 event . 
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If ICOT equals 0, SUBE3 calls TIMMY to determine the 
intercept position, single shot kill probability, SAM/SSM 
flight time, and the intercept possibility given the BLUE 
missile flight and warhead characteristics. If TIMMY re- 
turns the possibility flag -1, the intercept is not possible, 
but may be in the future. SUBE3 stores a Type 3 event for 
present time plus six seconds. If the indicator is zero, 
the intercept is not and will not be possible, and the sub- 
routine exits. If the indicator is one, the intercept is 
now possible and the program continues. If the intercept 
was possible, but the intercept point was in the hazard zone 
of another BLUE unit, TIMMY returns the indicator two. In 
this case, SUBE3 again stores a Type 3' event for present 
time plus six seconds. 

If the SAM/SSM requires the director for only the ter- 
minal phase of flight, SUBE3 determines whether the director 
can be scheduled for the intercept. If not, the subroutine 
stores a Type 3 event for the time the director will be 
free, minus flight time, plus terminal homing time, plus 
0.06 seconds. 

If scheduling is not needed,- or the director will be 
available when needed the missile is considered launched. 

The SAM-ASM- RAIDER and SAM tables are filled in and missiles 
on board and in the fleet counters are decremented. If 
there is coordination, and the firing ship is coordination 
capable, the target availability is up-dated. The launcher 
reload time is determined and stored. The channel 
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availability indicator is set to 1 and a Type 5 event is 
stored for intercept time. 

Prior to exit, if X(9) is set, SUBE3 prints out the 
firing ship number, the target number, and the reason for 
exiting the subroutine are printed. The breakdown for the 
REASON code for SUBE3 follows: 

0 Blue Salvo fired. 

1 Firing ship sunk or out of action. 

2 Target being engaged. 

3 Target already destroyed. 

4 Firing ship out of missiles. 

5 Director required for entire missile flight 
and is not available. 

6 Firing ship already engaging the target. 

7 Missile requires director only part time but 
all channels are busy. 

8 NONO table prohibits engagement. 

9 Launcher not available. 

10 ASM will impact prior to SAM launch. 

11 Director cannot be scheduled for the intercept. 

12 Director operation in horizon scan mode, busy. 

13 Target will not close BLUE missile envelope or 
TIMMY indicates no intercept. 

14 TIMMY indicates target may enter envelope later. 

15 TIMMY indicates target intercept point is too 
close to another BLUE unit. 

11. SUBROUTINE SUBE4 

SUBE4 determines whether the RED missile failed in 
flight, and if not, increments certain counters indicating 
ship hits and possible sinking. SUBE4 is called by EFFECT 
and it in turn calls RANDOM, and, if X(55) is set, GETIMX 
on entry and exit. 
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SUBE4 first breaks down the input code into the RED 
ASM/SSM number, and the BLUE target ship number. The code 
enters SUBE4 as the integer RRRBB. 'It then checks to see 
whether the ASM/SSM has been shot down, or had its launch 
negated due to the loss of the launching raider prior to the 
actual missile launch time. If neither is the case, the 
probability of successful flight is drawn from the X array 
and RANDOM is called to produce a random number. If it is 
determined that a failure occurred, the ASM table is set to 
3 and the subroutine exits. Otherwise the target ship is 
assumed to have been hit and the appropriate counters are 
incremented. If X(5) is not set, the subroutine compares 
the number of hits taken to the amount of damage the ship 
can withstand, and if too man; 7 bits have been taken, de- 
clares the ship sunk. 

12. SUBROUTINE SUBE5 

SUBE5 evaluates the BLUE missile intercept and performs 
checks on other ASM/SSM' if the target shot down or sunk was 
a RED raider. SUBE5 is called by EFFECT and it in turn 
calls STOREV, REACT, RANDOM, and if X(55) is set, GETIMX on 
entry and exit. 

SUBE5 is called with an input variable of salvo number. 
Using this, SUBE5 gets the target number, firing ship number, 
missile system, director, and channel number from the SAM- 
ASM-RAID table. The director reassignment and reacquisition 
delays are determined and stored, fleet and individual ship 
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statistical data are taken, and a check is made to insure 
that the target is still engageable. If the target has al- 
ready been sunk or shot down, the 'SAM-ASM-RAID table is set 
to -1, indicating an overkill. If the target is an ASM/SSM 
that failed in flight or has already impacted, the salvo is 
determined to have been a failure and the SAM-ASM-RAID table 
is set to zero. If the target is engaged, the probability 
of BLUE missile successful flight and support is pulled 
from the X array, and the single shot kill probability is 
removed from the SAM-ASM-RAID table. RANDOM is called to 
provide a random number. Success is determined by comparing 
PK to the random number drawn. This is done by determining 
PK as follows: 

PK = PS x ( 1- ( 1-PKSS ) n ) 

Where PS is the probability of successful flight, PKSS 
is the single shot kill probability, and N is the number of 
BLUE missiles in the salvo. If the random number is less 
than PK, the shot was 'successful . If the shot was unsuccess- 
ful and the missile requires a director, the channel is re- 
leased and Type 3, shoot, event is stored for present game 
time. If the missile doesn't require a director, just the 
shoot event is stored, and the subroutine exits. 

If the shot was successful and the target was an ASM/ 
SSM, the SAM, SAM-ASM-RAID, ASM, and kill counters are up- 
dated, a Type 6, channel free, event is stored, and the 
subroutine exits. If the target destroyed was an aircraft, 
or if the number of hits has been sufficient to disable a 
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surface raider, the following steps are taken prior to the 
regular table update and exit. A check is made to determine 
whether the raider has any ASM's o*r SSM's aboard. If so, 
the ASM counter is incremented and the appropriate numbers 
are entered in the ASM table indicating that these missiles 
were lost with the launching platform. Next the ASM table 
is checked to locate any missiles launched by the raider 
just destroyed. When one is found, its launch time is com- 
pared to raider kill time, and if the launch time is greater, 
entries are changed in the ASM table to correspond. Upon 
completing the search, the program continues to exit where 
it checks X(55), X(56), and X(58) for debug output. 

13. SUBROUTINE SUBE6 

SUBEG releases the fire control radar channel at the 
time of the event. SUBE6 calls GETIMX upon entry and exit 
if X(55) is set. 

The input value to SUBE6 is an integer code AABCD, 
where AA is the firing ship number, B is the missile system, 

C is the director number, and D the channel number. After 
breaking down the code, SUBEG sets the appropriate portion 
of the CHANEL array to zero, and exits. 

14. SUBROUTINE TIMMY 

TIMMY computes the predicted intercept point and flight 
time for the BLUE SAM/SSM and determines whether the inter- 
cept lies within the intercept contour. TIMMY is called by 
SUBE3 and TIMMY in turn calls TMDASM, ASMTIM, RADIAN, ZLNKA, 
LOCTIM, and if X(55) is set, GETIMX on entry and exit. 
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Inputs to TIMMY are target number, firing ship number, 
missile system number, present target position, and present 
game time. Outputs are single shot kill probability, flight 
time, and an .indicator flag which says the intercept is pos- 
sible, is not and will not be possible, or is not but may be 
possible in the future. 

Using the target number, TIMMY decides whether the tar- 
get is a raider or an ASM/SSM and determines the end of game 
position for raiders and target ship position for an ASM/SSM. 
If X(8) is set, firing sectors have been assigned and TIMMY 
checks to ensure that the assigned firing sector and the 
launcher train limits permit a launch. If launch is pro- 
hibited, REASON is set to 2, flag IND is set to zero, and 
the subrm.it. inp exits if the launch is possible, there ma.y 
be as many as two firing zones. These two occurrences are 
shown in Figure Bl. 
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Using the results of the ZLNKA reordering of the sector 
limits and launcher train limits, TIMMY determines whether 
there are two zones or one, and which of the bearings con- 
stitute the actual firing zone, or zones, limits. If no 
sectors have been assigned to the fleet, or this unit, the 
firing zone limits coincide with the launcher train limita.- 
tions . 

TIMMY next checks target altitude against the maximum 
SAM altitude. If the target is too high, REASON is set to 
eleven and the subroutine exits. Otherwise, the intercept 
parameters are determined. These intercept parameters are 
the values entered into the X-array that approximate the 
actual intercept contour. The intercept parameters are A, 

B, N, SLOPE, and SEEK. The actual intercept contour depends 
on target size, speed, and altitude as well as missile pa- 
rameters. There are usually more than one contour for any 
given missile system. The actual contours are found in the 
various technical publications dealing with the specific 
missile. OP-3594, Performance and Installation Character- 
istics Surface Missile Weapon Systems, also contains general 
contour nomograms, and the Applied Physics Laboratory of 
Johns Hopkins University has developed computer simulations 
which generate contours. A typical contour, along with the 
TIMMY approximation, is shown in Figure B2 . A is the off- 
range capability, B is the down-range capability, and N is a 
degree of curvature. 
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Figure B2a 

The Actual Intercept Contour 




Figure B2b 

The Approximation of the Contour 
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The maximum missile range for the given target altitude 
is determined, and if the target will not enter this range, 
as determined by TMDASM, TIMMY exits after setting REASON to 
seven. If the target is an ASM/SSM and the target ship is 
the firing ship, target course is found and SAM flight time 
is determined by: 

FT = DLT/ (SSPD+VEL) 

where DLT is the distance between the firing ship and the 
target, SSPD is the speed of the SAM, and YEL is the speed 
of the target. In this case the PK that is returned to 
SUBE3 is that for a non-crossing target. 

If the target ship is not the firing ship, or the tar- 
get is a raider, the following angles are determined: ALPHA, 
the angle between the target flight path and the line of 
sight at launch, and THETA, the angle in Cartesian coordi- 
nates of the present line of sight. The quantity GSINE is 
determined by: 

GSINE = SIN (ALPHA)* (VEL/SSPD) 

If the absolute value of GSINE is greater than one, the 
intercept is impossible and TIMMY exits after setting REASON 
equal to ten. Otherwise, the angles GAMMA, BETA, and DELTA 
are computed. 

GAMMA = ARCS IN (GSINE) 

BETA = 180° - GAMMA - ALPHA 
DELTA = 180° - BETA 
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GAMMA is the angle between the SAM/SSM flight path and the 
line of sight, BETA is the angle between the SAM/SSM flight 
path and the target flight path, and DELTA is the compliment 
of BETA. If X(59) is set, the values for target number, 
firing ship number, ALPHA, BETA, GAMMA, present target posi- 
tion, GSINE, and present range to the target are printed. 

PIC for this type of target is determined by DELTA. If DELTA 
is greater than 45° , the target is declared a crossing tar- 
get and the appropriate PK is retrieved from the X-array. 
Otherwise, the probability of kill for a non-crossing target 
is retrieved. Flight time is determined by: 

FT = DET * SIN (ALPHA ) / ( SSPD * SIN (BETA) ) 

Using present time plus flight time for the input, TIMMY now 
calls LOCTIM or ASMTIM, as appropriate, to determine the 
intercept coordinates. The range to the intercept is com- 
puted and compared to the minimum missile range. If the 
range is too short, IND is set to minus one, REASON to four, 
and the subroutine exits. The launch bearing is calculated 
and the first of two possible sets of firing zone bearings 
are compared to determine whether the launch bearing is with- 
in the firing bearings. If it is not, TIMMY checks to see 
if there was a second firing zone, if so it checks this one, 
otherwise, IND is set to minus one, REASON to three, and 
TIMMY exits. If the intercept angle is good, the intercept 
parameter is checked. A, B, and N are combined to form the 
equation : 

N N 

(X/A) + ( Y / B ) = 1 
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A value of two for N results in a circle or an ellipse. A 
value of* six results in a near rectangle. Fractional values 
for N are permitted. 

The "e3.1ipse" is translated so that its center is on 
the firing ship, and rotated so that the positive Y axis is 
opposite to the line of flight of the target. The intercept 
coordinates are substituted into the equation, in absolute 
value, and if the resulting number is less than or equal to 
one, the intercept will be within the ellipse. The next 
check is the seeker limit test. The y intercept value and 
the absolute value of the x coordinate are entered into the 
equation : 

Z = YI - SLOPE * XI 

If Z is greater than or equal to SEEK, the problem continues. 
If not, TIMMY exits after setting REASON to six. The sub- 
routine now checks the X-vector to determine whether there 
are any restrictions on how close an intercept can be to 
another friendly unit. If there are, the hazard zone is 
constructed. If the intercept falls within the hazard zone, 
REASON is set to nine, IND is set to two, and TIMMY exits. 

If the problem has progressed this far, the intercept 
is considered possible and the BLUE missile salvo is launched 
while IND is set to one. If the target is an ASM/SSM and 
the firing ship is the target, the program exits. If not, 
the target altitude is figured into the flight time prior to 
return to SUBE3. If X(9) is set, TIMMY prints out the ship 
number, target number, IND value, and REASON prior to exit. 
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Reason codes a.s used in TIMMY 



'1 Not used. 

2 No firing zones for this system. 

3 Firing bearing outside sector. 

4 Target too close. 

5 Target intercept point outside engagement contour. 

6 Target beyond seeker limits. 

7 Target beyond range and outbound. 

8 Intercept after ASM/SSM impact. 

9 Target too close to BLUE unit at intercept. 

10 GSINE greater than one. 

11 Target too high. 

15. SUBROUTINE JAMM 

The subroutine JAMM determines th-e range at which a 
monostatic radar can detect a given target in a multiple 
standoff and self-screening jamming environment. 

The calling routine, either SUBE1 in the case of raiders, 
or SUBE2 in the case of missiles, first determines that 
there are jammers in the problem and that at least one jam- 
mer is on the same frequency band as the detecting radar. 

If both situations obtain, subroutine JAMM is called. 

JAMM initially checks to see if the target will close 
the detecting unit within the minimum of the clear environ- 
ment, radar horizon, and theoretical range of the radar. If 
the target will not close within this minimum range, the 
subroutine is exited with a flag set to show that detection 
would not occur. In the case where the target will enter 
this minimum range, JAMM computes the time at which the tar- 
get will arrive at the closest point of approach (CPA) to 
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the detecting radar. Appropriate checks are made for sta- 
tionary targets and targets with only minor changes in their 
X or Y position components. JAMM next obtains the bearing 
from the radar to the target at the CPA position for use in 
selecting the appropriate side-lobe coefficient. Individual 
radars are then checked to see if they have been shot down 
prior to the target arriving at CPA* (applicable only for 
missile detections), if they are jamming on the frequency 
band of the detecting radar, and if they are above the de- 
tecting radar's horizon. If all of the conditions are met, 
the bearing from the jammer to the jammer's jamming target 
is found in order to determine whether the detecting radar 
is within the main beam or the side-lobe of the jammer. 
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jammer is then computed, taking into consideration the power 
of the jammer, the range from the jammer to the detecting 
radar, and the side-lobe/main beam relationship between the 
detecting radar and the individual jammer. The jamming power 
arriving at the detecting radar is then summed over all of 
the jammers and used to compute the burn-through range for 
the radar against the specified target. Prior to exiting 
the subroutine, the range to the target's CPA to the detec- 
ting radar is checked to ensure that it is less than the . 
burn-through range. If the distance to CPA is greater than 
the burn-through range a flag is set indicating that detec- 
tion did not occur due to CPA considerations. Otherwise the 
subroutine returns as the detection range the minimum of the 
clear, jamming, radar horizon, and theoretical ranges. 
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16. SUBROUTINE BIDET 



This subroutine computes the number of detection zones, 
ranges, and times of detection for a bistatic radar in 
either a clear or jamming environment. BIDET is called by 
SUBE1 and SUBE2, and it in turn calls GET I MX on entry and 
exit if X(55) is set. 

BIDET begins by setting the initial values of aircraft/ 
ASM start position, final position, altitude, and velocity. 
The transmitter and receiver positions and radar antenna 
height are established for the bistatic system under consid- 
eration. A check for non-moving targets is conducted prior 
to entering the closest point of approach (CPA) routine. In 
the case where the target has moved less than 0.1 nautical 
nii] 0 in both t h 0 X and Y diroction, t ho targot’s speed is 
considered to be zero and the target enters a special loop 
to determine whether it is detectable at its initial posi- 
tion within the constraints of both radar horizon and radar 
signal strength. 

For the moving target, the CPA of the target’s projected 
track to the transmitter is determined by solving simultane- 
ously the equations of the target's track and the line per- 
pendicular to the target's track and passing through the 
transmitter location. The CPA to the receiver is solved in 
a similar manner. If movement in either the X or Y direc- 
tions is less than 0.1 nautical mile, the CPA is established 
by setting the CPA coordinate for the small changing axis to 
the target's original position corresponding coordinate and 
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the alternate coordinate to the corresponding alternate co- 
ordinate' of the transmitter or receiver. For example, if 
the X coordinate of the target's positions changes less than 
0.1 nautical mile, the CPA to the receiver is set to the X 
coordinate of the target's initial position and the Y co- • 
ordinate of the receiver's position. 

The radar horizon for the transmitter and the receiver 
are then compared to the distance to their respective target 
track CPAs to determine whether the target enters the radar 
horizon of both the transmitter and the receiver. If it 
fails to enter both, the routine is exited with a no detec- 
tion. Otherwise, the times that the target first appears 
above and disappears below the transmitter and receiver ra- 
ri pr* tioi'izons sii * 0 computed. The time intepvsil th?*t the t si u — 
get is simultaneously above both horizons is determined by 
taking the maximum of the appearance times and the minimum 
of the disappearance times. If the target is never simul- 
taneously above both horizons; or the appearance time is 
greater than maximum game time; or the disappearance time is 
prior to start of game time, the subroutine is exited with a 
no detection. 

Next, detection times for the target based on radar 
signal strength are computed, and the intersection of these 
time intervals with the above joint radar horizon time in- 
terval determines the time during which the target is de- 
tected. No detection Occurs if the intersection of the two 
time sets (i.e. radar strength detection and radar horizon 
detection) is the null set. 
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If a detection has occurred in the clear environment, 
the subroutine checks to see whether any jammers are in use 
and if they are on the same frequency band as the bistatic 
radar. If jammers are in use, the time that the target 
arrives at its midpoint in the bistatic detection envelope 
is determined. This time is used in arriving at an average 
jammer position to compute the jamming power at the radar 
during the detection phase. This will tend to provide a 
somewhat conservative estimate of detection ranges in the 
face of jamming because the jammers are in general going to 
be nearer to the radar and the range square losses of the 
jamming power will be reduced. In the case of the station- 
ary target, midpoint time is set to an arbitrary three min- 
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start of game time to 
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the jammers toward the force. The subroutine then adds the 
jamming power from each jammer. It does this by checking 
all of the raiders to determine if the individual raider has 
been shot down prior to the midpoint time, and if the jammer 
is above the radar horizon and is jamming on the frequency 
band of the radar. If these conditions are met, the sub- 
routine examines the geometry of the jammer's position and 
the radar to determine the relationship between the jamming 
and radar antennas. With this relationship, the appropriate 
side- lobe/main beam antenna coefficients are used to adjust 
the jamming level. 

The constant in the bistatic range equation is then ad- 
justed to reflect the jamming environment and the computation 
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previously conducted in the clear environment are repeated 
to determine whether detection will occur in the jamming 
situation . 

The subroutine exits with the number of detection zones 
which occurred; the time and range of initial detection in' 
the first zone, the time of last detection in the first 
zone; and, if applicable, the time and range of first detec- 
tion in the second zone along with the time contact is lost 
in the second zone. 

17. SUBROUTINE DTI HD 

Subroutine DTIMD is a double precision subroutine that 
determines the real and imaginary roots of the fourth order 
polynomial which is used to solve for the bistatic radar de- 
tection times. Trie subroutine employs the Newton-Raphson 
iterative technique with the final iteration on each root 
performed using the original polynomial rather than the re- 
duced polynomial to avoid accumulated errors in the reduced 
polynomial. The subroutine was obtained from the Naval 
Postgraduate School Computer Center library. 

18. SUBROUTINE CPA 

Subroutine CPA obtains the closest point of approach of 
a target's track to a reference point by solving the simul- 
taneous equations of the line representing the target's 
track and the line through the reference point and perpen- 
dicular to the target's track. Checks are conducted for 
target tracks which pass through the reference point, 
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represent stationary targets, or targets which have only 
minor changes in either the X or the Y direction. 

19. SUBROUTINE REACT 

Subroutine REACT allows the use of reaction time delays 
through the use of probability distribution functions. This 
subroutine is similar to, if not identical to, one developed 
for use in FLOATS by APL at the Johns-Hopkins University. 

The distribution functions may be entered for seven classes 
of delay. They are: 

1 Target detectability to actual detection. 

2 Detection to ID and designation to the weapons 
direction equipment. 

3 Director assignment delay. 

4 Director acquisition delay. 

5 Firing delay. 

6 Target detectability to detection for horizon 

scan mode. 

7 Decision delay in horizon scan mode. 

Within these class areas, up to five separate types may 
be designated. REACT uses the type and class numbers to de- 
termine a delay time which it returns to the calling sub- 
routine. REACT is called by SUBE1 , SUBE2 , and SUBE5 . 

REACT in turn calls only RANDOM. 

The inputs to REACT are the class number and type, and 
the output is the delay in minutes. The class number is. 
used to locate that portion of the X-array where the delay 
is stored. The type number is used to localize the address 
within the area. In employing its Monte Carlo technique, 
REACT calls RANDOM and converts the returned uniform (0,1) 
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number to an integer. If this integer is 10, the subroutine 
sets it equal to 9. If the integer is zero, the RANGE, de- 
fined as the difference in time between two adjacent blocks, 
is set equal to the value in the first block. Otherwise, 
the address is determined using the class and type numbers 
and then adding on the integer. The delay consists of two 
parts. The first part is the time in the address block. 

The second part is determined by drawing another random num- 
ber and multiplying it by RANGE. Except in the case where 
the integer is zero, RANGE is determined by subtracting the 
time in the address block from that in the next highest 
block. In the "zero" case, the total delay consists of this 
fraction of RANGE. In all other cases, the fraction of 
RANGE is added to the time stored in the address block to 
obtain the total delay. At this point REACT returns to the 
calling routine. 

As an example of a distribution function, consider the 
waiting time for an aircraft for a catapult launch. Data is 
available for one hundred aircraft. The maximum waiting 
time was seven minutes. The distribution of waiting times 
was as follows: 

length of wait 0-1 1-2 2-3 3-4 4-5 5-6 6-7 ~ 

frequency 10 40 15 15 5 5 10 

Figure B3 

An Example of Waiting Time 
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The plot of length of time waiting against number of occur- 
rences would look like this: 




Figure B4 



Histogram of Waiting Times 



The probability of waiting between 0.0 and 1.0 minutes is 
0.1, the probability of waiting between 1.0 and 2.0 minutes 
is 0.4. When the histogram is transformed into a cumulative 
probability distribution function, the following results. 




Figure B5 

Cumulative Distribution Function of the Waiting Times 
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This distribution is what is entered in the X-array. 

The first block would contain the time for .1 cumulative 
probability. This would be one minute. Block two that for 
.2 cumulative probability, or 1.25 minutes, and so on up to 
block ten which would have entered a seven for seven minutes. 

20. SUBROUTINE PRIOR 

Subroutine PRIOR is called by SUBE3 . PRIOR calls REMOVE, 
STOREV, LOCTIM, ASMTIM, TARGET, and if X(55) is set, GET I MX 
on entry and exit. The purpose of PRIOR is to determine the 
most threatening target among those that a system may choose 
from at present game time. 

PRIOR first computes the value CHECK for the target 
being compared. CHECK is defined as maximum missile system 
range divided by present target range. TARGET is called to 
provide PRIOR with a candidate for comparison. The same 
measure of threat is obtained for the candidate ta.rget and 
stored in CHERRY. If the present value in CHECK is greater 
than or equal to CHERRY 1 , PRIOR goes to the next candidate, 
if there is one. If CHERRY is greater than CHECK, the val- 
ues of target number, channel number, and target position 
are saved. CHECK is made equal to CHERRY’ and PRIOR seeks a 
new comparison. 

If the value returned by TARGET in MORE is not zero, 
there is another candidate to be evaluated. . If MORE is zero, 
there are no more comparisons to make and the final section 
of PRIOR is entered. The stored value of most threatening 
target number is compared with the target sent to PRIOR 
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originally from SUBE3. If they are the same, the subroutine 
exits. If PRIOR has found a more threatening target, REMOVE 
is called to take the type 3 event associated with the new 
target out of the events list. STOREV is then called to 
enter a type 3 event for the target being intercepted. 

The debug switch X(9) is checked at exit and if it is 
set, old target, new target, missile system (in the code 
AABC were AA is ship number, B is system number, and C is 
director number), old channel, and new channel are printed. 
X(55) is checked and PRIOR returns to SUBE3. 

21. SUBROUTINE RCHECK 

Subroutine RCHECK evaluates the prospective target ship 
to determine if the raider can in fact "see" the target. 
RCHECK is called by SUBE2 and RCHECK calls GET I MX on entry 
and exit if X(55) is set. The inputs to RCHECK are the 
raider number and proposed target ship number. The two di- 
mensional distance between the two units is computed. The 
antenna heights on the target ship are compared and the 
tallest one is used to determine the ship's radar horizon. 

If the ship has no radar, a height of 15 feet is entered. 

The combined radar horizons, maximum theoretical radar range 
of the raider's radar, and the actual range are compared. 

If the target is within both the radar horizon and the maxi- 
mum radar range, the go/nogo flag is set to zero. If either 
test fails, the flag is set to one and the subroutine exits. 
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22. SUBROUTINE TMDASM 



This subroutine determines if the target, or other unit 
in question, v;ill enter within a certain input range of a 
BLUE unit, and if so,. at what time will it enter and leave. 
TMDASM is called by SUBE1, SUBE2 , and TIMMY. TMDASM in turn 
calls GET I MX on entry and exit if X(55) is set. 

TMDASM begins by determining the range from the refer- 
ence ship to the target unit. If this range is less than or 
equal to the input range, the time of entry is considered 
to be the input time. If the range to the target unit is 
greater than the input range, TMDASM determines the slope of 
a line representing target unit movement. Using this and 
the present target unit position as a reference point, TMDASM 
loi'iiis me iiae equucxoii. This equation and the equation for 
a circle with radius equal to the input range are solved 
simultaneously to determine the points of intersection. 

This is done by substituting the line equation into the 
circle equation and solving the quadratic equation. If the 
roots are imaginary or equal, the entry is either not made 
or the length of stay is infinitely small. These cases are 
rejected as no solution. Upon finding two real roots, the 
subroutine determines which of the two points the target 
unit would reach first, and which one it would reach second. 
If the target unit is already further down its track than 
the exit point, no entry is indicated. Otherwise the time 
of leaving the circle is determined and, if the target unit 
was not already inside the circle, the time of entering. 
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23. SUBROUTINE TMDBOM 

This subroutine determines when the raider in question 
will reach a given range from the force center, position 
(0,0). TMDBOM is called by MAIN, and JULIE. TMDBOM in turn 
calls GET1MX on entry and exit if X(55) is set. 

TMDBOM proceeds in the same way as TMDASM except that 
the reference position, center of the range circle, is al- 
ways force center instead of a firing ship. TMDBOM only 
computes the time that the raider will cross a given range. 
It neglects the departure time. 

24. SUBROUTINE ZLNKA 

Subroutine ZLNKA is a copy of the subroutine SHSORT, 
mentioned below, as it existed on 1 February 1975. ZLNKA 
inputs a vector of- real numbers and an ordered vector of 
integers. ZLNKA reorders the real vector from smallest to 
largest and reorders the integer vector to the same order as 
the real vector. ZLNKA is called by TIMMY. 

25. SUBROUTINE TJSORT 

Subroutine TJSORT is a modification of the subroutine 
SHSORT which is available from the Naval Postgraduate School 
Computer library. TJSORT is called by SUBE2 , and BIDET. 

This subroutine inputs a vector of real numbers and rear-, 
ranges them from smallest to largest. 
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26. SUBROUTINE LOCTIM 



LOCTIM computes the raider position at time T. It is 
called by MAIN, SUBE1, SUBE2 , SUBE3 , JAMM, BIDET, TIMMY, and 
PRIOR. LOCTIM in turn calls GETIMX on entry and exit if 
X(55) is set. 

LOCTIM first determines raider velocity and the time of 
the raider position stored in the RAID table. The distance 
the raider travels between the present time and the stored 
time (DL) is computed. If the velocity is zero, present tar- 
get position is set to the initial position and the program 
exits. The distance between the stored position and the end 
of game position (D) is determined and the ratio DL/D is com- 
puted. The X,Y position at the time desired is determined 

l-\ TT • 

‘-'J • 

XP - XZ - AX * DT 

YP = YZ - AY * DT 

where 

XP is the X position at the time of interest 

YP is the Y position at the time of interest 

XZ is the stored X position 
YZ is the stored Y position 

AX is the distance between XZ and the end of game 
position 

AY is the distance between YZ and the end of game 
position 

DT is the ratio DL/D 
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If the input flag IND has been set to 5, the new time 
and position are stored in the RAID table, otherwise the 
subroutine exits. 

27. SUBROUTINE ASMTIM 

ASMTIM computes the ASM/SSM position at time T. It is 
called by SUBE3 , TIMMY, and PRIOR. ASMTIM calls GET I MX on 
entry and exit if X(55) is set. 

The logic of ASMTIM is identical to that of LOCTIM but 
the definition of variables has been changed. 

XZ becomes the ASM/SSM X launch coordinate 

YZ becomes the ASM/SSM Y launch coordinate 

AX becomes the distance between the launch position 
and the X coordinate of the target ship 

AY bocomss “tho dist i?.nc o bstv/sor* "tbo l?„ i ui r icb posiirxon 
and tire Y coordinate of the target ship 

The ASM/SSM position is now computed in the same manner 
as in LOCTIM and the subroutine exits. The positions are 
not saved by ASMTIM. 

28. SUBROUTINE RANDOM 

Subroutine RANDOM is a generator of Uniform (0.0, 1.0) 
random numbers. It was programmed by Gerald P. Learmouth, 
of the W.R. Church Computer Center, Naval Postgraduate 
School. The following is taken from the library description 
of the subroutine: 
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"The basic generator is a Lehmer congruential 
of the form: 

X(N) = A*X(N-1) MOD (P) 

where A = 7**5 = 16807 and P = (2**31)-1. The 
resulting X(N)'s are uniformly distributed random 
integers between 1 and P. The uniform REAL*4 num- 
bers are formed by right-shifting the integer X(N) 
and appending a proper floating point exponent to 
form a number between 0.0 and 1.0." 

29. SUBROUTINE IIISTG 

Subroutine HISTG is found in the subroutine library of 
the Naval Postgraduate School subroutine library. A copy of 
HISTG has been included at the end of the program. The sub- 
routine is used if the user wants a statistical display of 
the number of hits a certain specified BLUE unit will take. 
Also provided is a statistical display of the time interval 
over which the RED missiles arrive, and one of the hits per 
minute. A histogram of the selected parameter is provided in 
all cases. HISTG is called by MAIN just prior to run com- 
pletion if X(85) has been set and the conditions detailed in 
section V.B have been met. 

30. SUBROUTINE RAIDIS 

This subroutine provides a plot of the raider tracks 
from their initial position to launch position. The maximum 
size of the plot is 800 nrn from force center. In the case 
where a raider does not launch weapons, the 30 minute game 
time position is used as the terminal point of the track. 

If the launch or 30 minute position is more than 800 nm from 
force center, a label of "no launch" is displayed at the 
raider's initial position. 
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31. SUBROUTINE DISPL 

DISPL is similar to RAIDIS except that the BLUE unit 
positions are plotted and they must' be within 200 nm of 
force center. A 100 nm range circle is included in the dis- 
play. 

32. SUBROUTINE SETIMX 

This subroutine sets the internal clock to zero. The 
entry GETIMX gives the internal clock time when called. 

33. FUNCTION RADIAN 

The function RADIAN converts degrees to radians. It is 
called by MAIN, JULIE, and TIMMY. The conversion is done 
with : 

RADIAN = RD * 0.017453293 

where RADIAN is the output angle in radians and RD is the 
input angle in degrees. 
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APPENDIX C 



RADAR DETECTION RANGE INPUTS 

This appendix provides the equations used to determine 
the 3 and a inputs for the individual radar classes. 

1. CLEAR ENVIRONMENT 

3 is the range at which a one square meter target can be 
detected in a non-jamming environment. It is a constant of 
the radar system and in the monostatic case provides the 
range of detection for the radar by the relation: 

JL 

Range = 3(e)' 1 where a is the target's radar cross section. 

The general equation for 3 is: 

P t G 2 A 2 i 
3 = Kp-)( 47T ) 3 ] 
m 

where P^_ is the power of the transmitter in watts; G is the 
antenna gain coefficient; A is the radar wave length in 
meters; and P is the minimum detectable signal of the radar 
with appropriate consideration for probability of detection 
and propagation losses. If the above units are used in 
calculating 3, the result will be a detection range in meters 
which must be converted to nautical miles for use in the 
simulation. Values of 3, in nautical miles for a one square 
meter target, are reported in Ref. 9 for standard U.S. Navy 
radar systems. 
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In the case of the bistatic system, the equation for 8 



becomes 



8 = 



P 4 G , G A 2 
r ( — — ) A- L ... 

Kp m (4iD : 3 



x 

4 



where the possibility of different gains of the transmitter 

and receiver must be considered and the propagation losses 

in P must take into consideration the differences in the 
m 

transmitter to target path and the target to receiver path. 
Additionally, in the bistatic situation, the detection 
range is not only a function of the range from the receiver 
(R ) to the target but it is also a factor of the range from 
the transmitter (R^) to the target. Thus the bistatic range 
equation becomes 

R, R = 8 2 oP 
t r 



The monostatic radar cross section provides a suitable esti- 
mate for the bistatic radar cross section [Ref. 10, p. 592]. 



2. JAMMING ENVIRONMENT 

a is the range at which a one square meter target can be 
detected in a jamming level of one watt/mhz. The burn through 
range for a monostatic radar in a multiple-standoff jamming 
environment is : 

R bt = oA(o/SJP)* 

where a is the target's radar cross section and SJP is the 
sum of the jamming power arriving at the radar. The jamming 
power from an individual jammer on the radar's frequency is 
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a function of the power of the jammer; the side-lobe/main 
beam relationships between the jammer and the radar; and the 
square of the range from the jammer to the radar. The equa- 
tion is: 

SJP = E C.JP./R? 
i 1 1 1 

where Ch is the side-lobe coefficient for jammer i and the 
radar; JP^ is the jamming power of jammer i on the radar's 
frequency; and iL is the range from the jammer to the radar, 
a is found by the relation: 

P G,^ G * i 

, t t r s 

a ( ‘4TTB n (S/J)' 1 

where P is the power of the transmitter; G is the gain of 
the radar (In bistatio systems G is the square root of the 
product of the transmitter and receiver gains.); B n is the 
receiver bandwidth of the radar; and S/J is the signal to 
jamming ratio for the desired probability of detection and 
appropriate propagation losses. In the bistatic range equa- 
tion in a multiple standoff jamming environment, the burn 
through range is found in a manner similar to the clear en- 
vironment to be 

R bt R t = a(a/SJP)" 

where is the burn through range, is the range from 

the bistatic transmitter to the target and the remaining 
factors are as defined above. 
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APPENDIX D 



PROGRAM ARRAY LISTING 



This appendix provides an alphabetical listing of the 
principal arrays used in the computer program. The listing 
contains the values of the arrays for the indicated array 
arguments. The arguments for the arrays are specified in 
the heading of their respective listing. The dimensions of 
the individual arrays are stated in the COMMON statement at 
the beginning of each subroutine. 



ASMT (I,J) 



The ASMT array contains data concerning the RED missiles. 
I AMS/SSM launch sequence number 
J=1 ASM/SSM status code as follows: 

0 = Launch failure 

1 = Successful launch 

2 = Killed by SAM 

3 = Inflight failure 

4 = Hit target ship 

5 = Launch platform shot down or sunk 

prior to launch 

J=2 Time of impact on target ship 
J=3 Raid number of launch platform 



J=4 Target ship number 



J=5 Time of launch 



J=6 Weapon class 



J=7 X coordinate of launch position 



J=8 Y coordinate of launch position 
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CHANEL ( I , J , K, L) 


An array 


used to keep track of which channels are busy. 


I 


Ship number 


J 


System number 


K 


Director number 


L 


Channel number 

0 = Free 

1 = Busy 




JACK (I,J) 



This array stores the engagement and detection statistics 
for individual ships. 



I 


Ship number 


J=1 


Chances to engage. The number of targets that 
will enter within the maximum missile range of 
the ship 


J=2 


Number of engagements 


J=3 


Number of targets destroyed 


J=4 


The number of raiders that will enter the mis- 
sile envelope. 


J=5 


The number of raiders that will enter the en- 
velope and are detected 


J=6 


The number of ASM/SSM that will enter the 
envelope 


II 


The number of ASM/SSM that will enter the 
envelope and are detected 




NADAR (I,J) 



This array indicates the operational status of the on- 
board BLUE equipment. 

I Ship number 
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J=1 



Search radar number 1 



J=2 

J=3 

J=4 



J=5 

J=6 

J=7 

J-8 

J-9 

J=10 

J=ll 

J~12 

J=13 

J=14 

J=15 

J=16 



J=17 

J=18 

J=19 



0 = down 

1 = up 

Search radar number 2 
Search radar number 3 

Battery one, fire control radar number 

0 = down 

1 -• operational normal mode 

2 = no director needed for missile 

3 = operational in horizon scan mode 

Battery one, fire control radar 2 
Battery one, fire control radar 3 
Battery two, fire control radar 1 
Battery two, fire control radar 2 
Battery two, fire control radar 3 
Battery three, fire control radar 1 
Battery three, fire control radar 2 
Battery three, fire control radar 3 
Battery four, fire control radar 1 
Battery four, fire control radar 2 
Battery four, fire control radar 3 

Launcher for battery one 

0 = non operational 

1 = operational 

Launcher for battery two 
Launcher for battery three 
Launcher for battery four 



1 
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RAID ( I , J ) 



This array stores information on the raider's location 



and weapons. 




I 


Raider number 


J=1 


X vector starting index of the raid class 


J=2 


Time that raid was at position indicated in 
J=3 , 4 


J=3 


X coordinate at time T 


J=4 


Y coordinate at time T 


J=5 


X coordinate at end of game 


J=6 


Y coordinate at end of game 


J=7 


X coordinate at launch if the raider has weapons 
and can reach the launch position, otherwise 
the X coordinate at end of game. 


J=8 


Y coordinate at launch, same as J=7 


J=9 


Raider altitude in feet 


J=10 


Raider's radar horizon 


J=ll 


Weapon class of raider's longest range weapon 


J=12 


Total number of weapons assigned 


J=13 


Status 

0 = destroyed or sunk 

N = number of hits the raider may take until 
destroyed 


J=14 


X coordinate of raider at end of game or X 
coordinate of target ship for ASM/SSM launched 
by raider I 


J=15 


Y coordinate comparable to J=14 


J=16 


Time raider will reach second launch point if 
it is carrying two types of missiles 


J=17 


Time the raider was destroyed 
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SAMASM (I,J) 



This array contains information covering the BLUE missile 
launch, target, engagement outcome, and probability of kill. 



I 


BLUE salvo number 


J=1 


Target number, 1-90 are raiders, 91-191 are 
ASM/SSM 


J=2 


Ship number followed by system number 


J=3 


Director number followed by channel number 


J=4 


Time of launch 


J=5 


Engagement outcome. Negative number indicates 
an overkill, zero a failure, and a positive 
number a successful shot. 


J=6 


Time the channel will be free for reassignment 


J=7 


Single shot probability of kill 


The SHIP 


SKIP (I,J) 

array contains information concerning the BLUE 



units 



I 


Ship number 


J=1 


Ship class 


J=2 


X coordinate of ship position 


Ch 

II 

CO 


Y coordinate of ship position 


J=4 


Number of hits ship can receive until declared 
out of action. If 0, ship is out of action. 


J=5 


Ship class starting index 


J=6 


Missile status code 

0 = no missile capability 

1 = missile capability 
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SMAT (I,J) 

The SMAT array lists the results of the SAM firings. 



I 


SAM number 










J=1 


SAM class 


index 










J=2 


Firing ship number 










CO 

ii 

■"D 


Number of 


missiles 


; remaining in 


the battery 




firing the 


! SAM 










J=4 


Time fire 


control 


radar 


channel 


1 


free 


J=5 


Time fire 


control 


radar 


channel 


2 


free 


J=6 


Time fire 


control 


radar 


channel 


3 


free 


J=7 


Time fire 


control 


radar 


channel 


4 


free 


J=8 


Time fire 


control 


radar 


channel 


5 


free 


J=9 


Time fire 


control 


radar 


channel 


6 


free 


J=10 


Number of 


SAMs in 


the salvo 







STATS (I,J) 

The STATS array is used to compile statistics concerning 
the BLUE units. 

I Ship number 

J=1 Sum of the detection ranges on raiders 

J=2 Sum of the detection ranges on missiles 
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APPENDIX E 



GLOSSARY FOR PROGRAM VARIABLES 



1. MAIN PROGRAM 

ALFA REAL*8 variable used for conversion of degrees to 

radians . 

AMF Small real constant used for conversion from float- 

ing point variables to integer variables. 

ASMT The table containing RED ASM/SSM launch and target- 

ing information. 

BWJ Beam width of the jammer. 

BWR Beam width of the radar receiver. 

HISA Vector used to store data for use with subroutine 

HISTG. 

HISO Vector used to store data for use with subroutine 

HISTG. 

HISTG Subroutine in the NPS computer library that pro- 
vides a histogram display of an input vector. 

HISX Vector used to store data for use with subroutine 

HISTG. 

IBAT An array used to keep track of the number of mis- 

siles remaining in each missile battery of a BLUE 
unit . 

IFAG Integer constant containing the number of RED mis- 

siles engaged by the BLUE fleet. 

IFBG Integer constant containing the number of RED 

raiders engaged by the BLUE fleet. 

IHIT Vector containing the number of hits received by 

each BLUE unit. 

I IF The number of RED missiles that failed in flight. 

IKS The number of RED missiles shot down by BLUE SAMs. 

ILF The number of RED missiles that failed when launched. 
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ILS 


The number of RED missiles that were in flight at 
the end of the iteration. 


IMF 


The number of RED missiles destroyed when the 
launching aircraft or ship was destroyed. 


IMK 


The number of RED missiles destroyed by SAMs over 
all iterations. 


I MM 


The number of RED missiles destroyed prior to launch 
over all iterations. 


INSECT 


The original seed value used with the random number 
generator. Used to reproduce the random number 
chain if X(54) is set to 1. 


IRT 


The starting address in the X vector for the next 
five input quantities. 


ISAG 


An array containing the number of RED missiles 
first engaged by each BLUE ship. 


ISBG 


An array containing the number of RED raiders first 
engaged by each BLUE ship. 


ITS 


The number of RED missiles that impacted on their 

t 9 yrrof 
w o ^ • 


JACK 


The array containing the BLUE ship detection and 
engagement counters. 


KOUNT 


The iteration number. 


MARY 


The vector containing the data for the hit fre- 
quency histogram. 


NADAR 


The array that stores the BLUE radar and missile 
system availability for each iteration. 


NASM 


The number of ASM/SSMs. 


ND 


Index used to convert true bearings to cartesian 
bearings . 


NI 


The number of iterations to be performed. 


NMIF 


The number of missiles in the BLUE fleet. 


NR 


The number of armed and unarmed RED raiders. 


NRTGT 


The vector containing the RED units target list 
assignment . 
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NS 


The number of BLUE ships. 


NSAM 


The number of BLUE missile salvos fired. 


PLIST 


The array containing the target lists. 


RAID 


The table containing the RED raider position and 
weapon information. 


RAT 


REAL*8 variable used to convert degrees to radians. 


RXATA 


Dummy array used to read in data and fill the X 
vector . 


SAM ASM 


The array containing the BLUE missile RED target 
engagement information. 


SET I MX 


The subroutine that sets the internal clock if 
X(55) is set to 1. 


SHIP 


The table containing the BLUE ship position infor- 
mation . 


SMAT 


The table containing the BLUE missile battery en- 
gagement data. 


STATS 


HP Ho oproy tho RT TTT7 nn-i t rlof np-f-i on rorm-no 

i -X -* *-«. — C-V J *..* W — — X — < O 4.-1 (-1 XI a. 1/ M V-— C/ O J. J.X Uxi \-< • 


T 


Time. 



2. SUBROUTINE EFFECT 



IUNIT 


The code containing the units involved in playing 
the next event. 


NEXTEV 


The next event type. 


TIME 


The time the next event is to be played. 


TIMES 


The internal clock time used for entry and exit 
if X(55) is set. 



3. SUBROUTINE JULIE 



BNG 


A REALMS variable used to convert degrees to radians 


CHANEL 


The array that keeps track of which BLUE missile 
channels are busy. 


IA 


The X vector address of the fire control radar 
class number. 
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IBAT 


The array that keeps track of the number of BLUE 
missiles remaining in each battery. 


IC 


The radar class of the fire control radar associated 
with the number one director of the indicated BLUE 
missile battery. 


IFAG 


Integer constant containing the number of RED mis- 
siles engaged by the BLUE fleet. 


IFBG 


Integer constant containing the number of RED 
raiders engaged by the BLUE fleet. 


IFIRST 


Flag used for checking a target's priority. 


IFLASK 


Flag used to set the time by which all RED raiders 
must reach their launch point in order to launch 
their missiles. 


I FOG 


Flag to indicate a) horizon search mode, and 
b) the system has a launcher that is operable. 


I GAGE 


A vector that keeps track of which targets have 
been engaged by BLUE missiles. 


IHIT 


A vector that keeps track of how many hits each 
BLUE unit has taken. 


IJ 


The ship class starting index in the X vector. 


I OR 


The raid class starting index in the X vector. 


IP 


Missile system class number. 


IS 


The ship class starting index in the X vector. 


ISAG 


Integer constant to determine how many RED missiles 
the ship was the first to engage. 


ISBG 


Integer constant to determine how many RED raiders 
the ship was the first to engage. 


IT 


Search radar class number. 


ITEMS 


The number of events that have been retrieved from 
the events list. 


I VENTS 


The array containing the events list. 


JACK 


The array containing the ship detection and engage- 
ment statistics. 
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JAM 

KFOG 

LOVER 

MAMBAT 

MISLFT 

NADAR 

NASM 

NONO 

NR 

NS 

NSAM 

PROB 

RAID 

RELOAD 

RGE 

SHIP 

SMAT 

STATS 

TGTAV 



The vector that keeps track of the number of jammers 
operating on a given radar frequency band. 

Flag to indicate that a search radar is available. 

The number of events in the events list. 

The number of missile batteries available to the 
BLUE force. 

The array that keeps track of how many missiles a 
RED raider has left on board. 

The array that keeps track of the operational mode 
of the BLUE equipment . 

The number of RED missiles. 

The array that keeps track of engageability of a 
RED unit or missile with respect to each BLUE unit. 

The number of armed and unarmed RED raiders. 

The number of BLUE units. 

The number of BLUE missile salvos. 

The probability that a given piece of equipment is 
operational . 

The table with the raider location and weapon in- 
formation . 

The array that stores the time the launcher will 
be reloaded by after a missile firing. 

The launch range of the longest ranged missile 
carried by each raider. 

The table containing the ship position information. 

The table with the BLUE missile system launch in- 
formation . 

The table of BLUE detection ranges. 

The table that tells when a target will be avail- 
able for engagement by another BLUE unit. 
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4 . 



SUBROUTINE EVENT 



ICODE The participating units-event type code that is 
stored in the events list. 

ID The event type. 

IT The time of the event converted to the integer 

format used to store it in the events list. 

ITB Used to determine where in the events list the next 

event to be stored should go. 

ITEMS The .number of events retrieved. 

ITIM The time as stored in the events list. 

ITM The time as stored in the events list. 

IUNIT The units involved in playing the next event. 

IVENTS The events list. 

IVET The participating units as stored in the events 

list . 

LOVER Th o -n 1 1 ni v> nr of events stored ^n tlfio events list. 

NEXTEV The next event type. 

NUMBER The code for the participating units. 

T The time of the event. 

TIME The time of the event. 

5 . SUBROUTINE REMOVE 

IUNYT The code stored in the events list of the event to 
be removed from the events list. 

IVENTS The events list. 

LOVER The number of items in the events list. 

6 . SUBROUTINE TARGET 

IFIRST The flag that indicates that the counter has been 
zeroed. 

IT The game time in events list code. 
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I UN IT 



I UN IT 


The code of the event to be compared by PRIOR. 


I VENTS 


The events list. 


MORE 


The flag indicating that there is an event with 
the same playing time and of type three. 


NPGS 


The counter that is used as an index in searching 
the events list. 


TYME 


Present game time. 



7. SUBROUTINE SUBE1 



ALT 


The raider's altitude. 


AREA 


The fourth root of the raider's radar cross section 


UELR 


The detection delay. 


FCRH 


The radar horizon of the fire control radar. 


FRG 


The minimum of the clear environment detection 
range, the radar horizon, and the radar's theo- 
retical maximum detection range. 


HFCSR 


The radar horizon of the fire control radar. 


HORZ 


The combined radar horizons of the radar and the 
target . 


HRAG 


The combined radar horizons of the radar and the 
target . 


HRFC 


The combined radar horizons of the radar and the 
target . 


HSFR 


The minimum of the clear environment detection 
range, the radar horizons, and the maximum theo- 
retical detection range. 


HT 


The height of the fire control radar antenna. 


IA 


Used to check the class numbers of the various 
systems . 


IC 


Missile system class number. 


ICH 


Used to enter the NADAR table. 


ICL 


Fire control radar class number. 
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I CPA 

IFA 

IFC 

IFCR 

IF 

IR 

ISI 

JACK 

JAM 

JAMM 

JCPA 

JFOG 

KEY 

KFLAG 

LCPA 

MFOG 

MJL 

MTY 

NADAR 

NASA 

NASM 

NB 

ND 



Flag telling whether the target will come within 
the maximum detection range. 

Integer to store the radar number. 

Integer to store the battery number. 

Fire control radar class number. 

Radar frequency band. 

The raid number. 

The ship class starting index in the X vector. 

The table of engagements and detections. 

The table that keeps track of how many jammers 
there arc on each radar frequency band. 

The subroutine for determining the detection range 
of a monostatic radar in a jamming environment. 

The flag that tells whether or not the target will 
come within the detection range. 

The flag that tells whether or not this is the 
second bistatic detection zone. 

The type of probability distribution associated 
with this unit's time delays. 

Flag that says that this unit has a search radar 
that is operational. 

Flag that tells whether or not a target will close 
to within the detection range. 

Counter giving the number of second bistatic de- 
tection zones to be processed. 

Index for table lookup. 

Missile system class number. 

The table containing the equipment availability. 
The flag indicating that raider is an aircraft. 

The number of ASM/SSM launched or destroyed. 

The battery number. 

The number of bistatic detection zones. 
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NF 

NFOG 

NJ 

NONO 

NSR 

NSURF 

PROB 

RAGJ 

RAID 

RCPS 

REACT 

RFC.T 

RG 

RGE 

RGF 

RGFC 

RGI 

RGIA 

RGIB 

RGIC 

RGM 

RIP 

RRFC 



The director number. 

Flag telling whether or not there are any more 
bistatic detection zones to check. 

The number of jamming raiders. 

The table that contains the reason for non- 
engageability . 

The search radar index. 

Flag that says that this ship has a surface to 
surface missile system. 

Probability of acquisition. 

Detection range in a jamming environment. 

The table with raider location and weapon informa- 
tion . 

Range at CPA squared. 

Subroutine for determining reaction time delays. 

Fire control radar detection range in the jamming 
environment . 

The minimum of the clear environment detection 
range, the radar horizon range, and the theoreti- 
cal maximum detection range. 

Detection range. 

Jamming environment detection range. 

Detection range of the fire control radar. 
Detection range in a jamming environment. 

Storage for second bistatic detection range. 
Storage for second bistatic detection range. 
Storage for second bistatic detection range. 
Missile range. 

The maximum detection range, or present target 
position, whichever is smaller. 

Clear environment fire control radar detection 
range . 
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RRH 


Radar horizon of the search radar. 


R2 


The detection range in the second bistatic detec 
tion zone. 


SALT 


The radar horizon of the raider. 


SAMVEL 


The BLUE missile's speed. 


SHIP 


The table with the ship position information. 


SIGMA 


The raider radar cross section. 


SRAG 


Missile range. 


SRFC 


Fire control radar clear environment detection 
range . 


STATS 


Ship's detection range storage. 


SVEL 


Missile speed. 


TD 


Time raider will close to the given range. 


TDES 


The designation time delay. 


TDET 


r PV->/*N -t-nmo. r>"F — ' om 

a v ^ j- 


TD2 


The time of detection in the second bistatic de- 
tection zone. 


TFIR 


The last time to fire. 


TFIRE 


The firing time delay. 


TFR 


Time the target will close to fire control 
acquisition range. 


THR 


Theoretical radar maximum detection range. 


THSCAN 


Horizon scan detection delay. 


TI 


Time the target will close to the given range. 


TIME 


Present game time. 


TIML 


Time the target will go beyond the given range. 


TIMLA 


Storage for lose detection time. 


TIMLB 


Storage for lose detection time. 


TIMLC 


Storage for lose detection time. 
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TIML2 


Time of lose detection in the second bistatic de- 
tection zone. 


TLS 


Last time to fire. 


TRGE 


Time of detection. 


TSHOOT 


Time of launch. 


TSR 


Total reaction delay. 


TTA 


Storage for detection time in second detection zone 


TTB 


Storage for detection time in second detection zone 


TTC 


Storage for detection time in second detection zone 


TWDS 


Reaction time delay for weapons designation equip- 
ment . 


VEL 


Raider speed. 


XCP 


The X coordinate of the CPA. 


YCP 


The Y coordinate of the CPA. 



8 . SUBROUTINE SUBE2 



ALT 


Altitude of the RED ASM/SSM during flight in feet. 


ALTN2 


Altitude of the RED ASM/SSM squared in miles. 


AREA 


The fourth root of the radar cross sectional area 
of the RED missile. 


ASMT 


The table containing RED missile launch and target- 
ing information. 


DELR 


The detection time delay. 


DELTAX 


The distance in the X direction between the target 
ship and the ASM/SSM launch point. 


DELTAY 


The distance in the Y direction between the target 
ship and the ASM/SSM launch point. 


DIST 


The three dimensional distance from the ASM/SSM 
launch point to the target ship. This distance 
assumes launch at cruising altitude. 


HORZ 


The combined radar horizons of the firing ship and 
the RED missile. 
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HR AG 


The combined radar horizons of. the firing ship and 
the RED missile. 


HT 


The BLUE radar antenna height . 


TA 


Fire control radar class number. 


IC 


The BLUE missile system class number. 


ICL 


Fire control radar class number. 


ICOP 


Constant for X vector entry set by SUBE2. 


I CPA 


Flag telling whether or not- a RED missile will come 
within the detection range. 


I DP 


Index counter used to make up the RED raider's 
temporary target list. 


I FA 


The director number. 


IFC 


The battery number. 


IFCR 


BLUE missile system class number. 


IFLAG 


Flag indicating RED unit is firing a two missile 

S3 1 VO . 


IFLASK 


Flag indicating the first air launch has occurred, 
set TMAX. 


IHIT 


Table of hits taken by BLUE ships. 


IJ 


Raider class number. 


IJI 


Radar frequency band. 


IOU 


Target ship number. 


IPQ 


Target hit number. 


IR 


Raider number. 


ISI 


Ship class starting index in the X vector. 


JACK 


Ship detection and engagement table. 


JAM 


Table that keeps track of the number of jammers on 
the radar frequency bands. 


JAMM 


Subroutine to determine monostatic radar detection 
ranges in a jamming environment. 
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JFLAG 


Flag indicating that there are more second bistatic ‘ 
detection zones to evaluate. 


JFOG 


Flag indicating that there were two bistatic detec- 
tion zones. 


KATHY 


Flag indicating that a second type ASM/SSM is being 
processed. 


KEY 


Type of probability distribution to be used in 
computing the reaction time delay. 


KFLAG 


Flag indicating BLUE unit has an operational search 
radar . 


KING 


Flag for one or two missile salvos. 


KQED 


Flag indicating whether or not raid could target a 
particular BLUE unit. 


LNUM 


The class number of the raider's long range missiles 


MARIE 


Number of hits required to put the BLUE unit out of 
action . 


MFOG 


Counter indicating the number of second bistatic 

2OD.0S . 


MISLFT 


Array to indicate the number of ASM/SSM left on 
board a raider. 


NADAR 


The array that stores BLUE equipment operability. 


NASM 


The number of ASM/SSMs processed. 


ND 


The number of bistatic detection zones. 


NEXT 


Flag indicating the short range RED missiles. 


NFOG 


Flag indicating there are more detection zones to 
process . 


NJ 


The number of jammers. 


NMIS 


The number of ASM/SSM of this type to be launched. 


NONE 


The index used to draw random numbers from the 
random number vector . 


NONO 


The array specifying the target-firing ship rela- 
tionship . 


NOW 


The type missile being processed now. 



NR 

NRTGT 

NS 

NTOT 

PLIST 

PROB 

QRT 

RAGJ 

RAID 

RCPS 

RELI 

RFCJ 

RG 

RGE 

RGF 

RGFC 

RGI 

RGIA 

RGIB 

RGIC 

RGM 

RGS 

RIP 

RRFC 



The number of raiders. 

'The target list assignments. 

The number of BLUE units. 

The total number of ASM/SSM on board. 

The target lists. 

The probability of acquisition in horizon scan 
inode of operation. 

The vector of (0,1) random numbers used to deter- 
mine missile launch times. 

The search radar detection range in a jamming en- 
vironment . 

The table with the raider location and weapon in- 
formation . 

The slant range to CPA. 

The RED missile's inflight reliability. 

r)of opHon range of fire control a cl a r i™ a jamming 
environment . 

The minimum of the clear, jamming, and theoretical 
detection ranges. 

Detection range. 

Detection range. 

Fire control radar detection range. 

Detection range. 

Storage for secondary bistatic detection range. 
Storage for secondary bistatic detection range. 
Storage for secondary bistatic detection range. 
BLUE missile range. 

Detection range squared. 

Minimum of present range and detection range. 
Horizon search clear environment detection range. 
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RRH 


Radar horizon. 


RT 


‘Cumulative probability of targeting this target or 
one previous. 


R2 


Second bistatic detection range. 


SALT 


The radar horizon of the RED weapon. 


SAMVEL 


BLUE missile speed. 


SHIP 


The table containing ship position information. 


SIGMA 


The RED weapon radar cross sectional area. 


SPD 


The RED weapon speed. 


SRAG 


Search radar detection range in a clear environment 


SRFC 


Fire control radar detection range in a clear en- 
vironment . 


SRGE 


BLUE ship missile range. 


STATS 


Array containing the sum of a ship's detection 
ranges. 


SVEL 


BLUE miss„ile speed. 


TAC 


Acquisition time delay. 


TD 


Time of detection. 


TDES 


Designation time delay. 


TDET 


Time of detection. 


TD2 


Time of detection in the second bistatic detection 
zone . 


TFIR 


Last time to fire based on search radar range and 
maximum range missile. 


TFIRE 


Firing time delay. 


TFR 


Time target will come within the fire control de- 
tection range. 


THR 


Maximum theoretical radar range. 


TUSCAN 


Reaction time delay for horizon scan mode of opera- 
tion . 
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TI 


Time target will come within the search radar de- 
tection range. 


TIFA 


Time target will come within the horizon scan 
detection range. 


TIME 


Present game time. 


TIMF 


Time target will exit the fire control radar detec- 
tion range. 


TIML 


Time the target will go beyond the search radar 
detection range. 


TIMLA 


Storage for second bistatic detection zone loss of 
target time. 


TIMLB 


Storage for second bistatic detection zone loss of 
target time. 


TIMLC 


Storage for second bistatic detection zone loss of 
target time. 


TIML2 


Time target will go beyond detection range in the 
second bistatic detection zone. 


TINT 


The time interval over which al ] missiles must be 
launched . 


TLS 


Last time to fire based on fire control radar range 


TM 


Proposed RED missile launch time. 


TMAX 


Time limit on launch of RED missiles. 


TRGE 


a) Hazard zone range squared, b) time of . detection . 


TSHOOT 


Time of proposed BLUE missile launch. 


TSR 


Total reaction time delay. 


TTA 


Storage for second detection zone detection time. 


TTARG 


The temporary target list. 


TTB 


Storage for second detection zone detection time. 


TTC 


Storage for second detection zone detection time. 


TWDS 


Delay for weapons direction equipment processing. 


TX 


Time between RED missile launchings. 


VEL 


RED missile speed. 
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XCPA 



XCPA 

YCPA 


CPA position in the X direction. 

/ 

CPA position in the Y direction. 


z . 


(0,1) random number. 



9. SUBROUTINE SUBE3 



ASMT 


The table containing the RED missile launch and 
targeting information. 


CHANEL 


The array that indicates which system channels are 
in use. 


EBTJR 


The Y coordinate of CPA. 


FT 


BLUE missile flight time. 


I BAT 


Array giving the number of BLUE missiles available 
by battery. 


I COT 


Flag indicating channel assigned. 


IP 


Ship class starting index in the X vector. 


IU 


Input code giving target, ship, system, and direc- 
tor information. 


JAPL 


A channel number. 


JULIA 


RED missile number for table lookup. 


K 


Firing ship number. 


KA 


Target ship number. 


M 


Number of original target. 


MFOG 


Flag used to check director scheduling. 


Ml 


Target number after searching PRIOR. 


NADAR 


Array to check launcher operability. 


NB 


Director number. 


NC 


Channel number. 


NCH 


Number of channels associated with this radar. 


NCI 


Channel number. 


NMS 


Missile system number. 
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NSAM 


BLUE salvo number. 


PROB 


Probability of kill. 


RAID 


Table of raider positions. 


REASON 


Code for non-engageability . 


RELOAD 


Table that tells if the launcher is loaded. 


RGE 


Systems maximum range. 


RJL 


The X coordinate of CPA. 


RLT 


Launcher reload time. 


SAMASM 


The table of BLUE missile and RED target engagement 
data. 


SHIP 


Table of ship positions. 


SMAT 


Table of BLUE missile system data. 


TCF 


Dummy vector for storing the times the channels 
will be free for other use. 


TFT 


Time director first available for channel use. 


TF 


Time director needed. 


TGTAV 


Time target available for engagement by another 
unit. 


TIM 


Time channel free. 


TIME 


Present game time. 


TL 


Time director no longer needed. 


TT 


Time director needed for homing. 


XP 


Present target X coordinate. 


YP 


Present target Y coordinate. 



10. SUBROUTINE SUBE4 



ID 


Used as a class indicator, i.e. ASM type, or as an 
integer constant, i.e. number of hits to sink the 
target ship. 


IHIT 


The table containing the number of hits sustained 
thus fax - in the game. 
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IU 


The input code containing the RED missile number 
and the BLUE target ship number. 


K 


The target ship number. 


NAM 


The RED ASM/SSM number for table lookup. 


NAS 


The RED ASM/SSM number in the input code. 


PB 


The probability of successful flight for the RED 
missile. Also used in determining the input and 
output time of the subroutine if X(55) is set. 



11. SUBROUTINE SUBE5 



ASMT 


The table containing the data on the RED missile 
engagement . 


CHANEL 


The table that keeps track of which BLUE missile 
channels are being used. 


ID 


An integer constant used for the starting address 
in the X-array of the BLUE firing ship's class 
inputs. Also used as the number of RED ASM/SSMs 
launched up to the present game time. 


I FAG 


Integer constant that totals up the number of RED 
missiles engaged by the BLUE fleet. 


IFBG 


Integer constant that totals up the number of RED 
raiders engaged by the BLUE fleet. 


IGAGE 


Integer constant that keeps track of which targets 
have been engaged. 


IP 


Radar class of the BLUE fire control radar asso- 
ciated with the engagement. 


IQ 


BLUE missile system class number. 


ISAG 


Integer constant that totals up the number of RED 
missiles first engaged by this firing ship. 


ISBG 


Integer constant that totals up the number of RED 
raiders first engaged by this firing ship. 


JACK 


The table that keeps track of the number of missiles 
and raiders engaged by the firing ship. 


JAM 


The table that totals up the number of jammers 
operating on a particular radar frequency band. 
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JFK 


Integer used in table lookup. 


K 


The BLUE firing ship. 


KEY 


As used in the statement -- CALL REACT (KEY,4,TA) - 
where 4 is the type of delay, KEY is the class of 
the particular distribution within the type, and 
TA is the delay. 


M 


The BLUE salvo number. 


MISLFT 


The table that keeps track of the number of mis- 
siles remaining on board the RED raider. 


N 


The number of missiles in the BLUE salvo. 


NA 


The RED ASM/SSM number for table lookup. 


NASM 


The number of RED ASM/SSMs processed thus far. 


NB 


The BLUE ship's director involved in the engagement 


NC 


The channel number of the fire control system en- 
gaged in the intercept. 


NDU 


The starting address of the class associated with 
the raider involved in the intercept. 


NMS 


The number of the ship's missile system involved 
in the intercept. 


NT 


The target number, 1-90 for a raider, 91-190 for a 
RED missile. 


NTOT 


The number of missiles left on board the raider. 


PB 


The probability that the BLUE missile will receive 
the proper support during flight. 


PROB 


The overall probability of a successful intercept. 


SAMASM 


The table that keeps BLUE salvo information. 


SMAT 


The table that keeps BLUE system information. 


TA 


The time delay for target acquisition. 


TF 


The time delay for firing. 


TK 


The time delay for kill assessment. 
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12. SUBROUTINE SUBE6 



CHANEL 


The table that keeps track of which BLUE missile 
channels are being used. 


K 


The firing ship. 


M 


The input code containing the ship, system, direc- 
tor, and channel number to be released. 


NB 


The director number. 


NC 


The channel number. 


NMS 


The system number. 


T 


The input and output time for the subroutine, used 
if X(55) is set. 



13. SUBROUTINE TIMMY 



A 


The cross range factor in the intercept contour. 


ALPHA 


The angle between the target flight path and the 
line of sight to the target. 


ALT 


The target altitude. 


ASMT 


The table containing the data on the RED missile 
engagement . 


ASMTIM 


Subroutine to determine RED missile position at a 
given time. 


B 


The downrange factor in the intercept contour. 


BETA 


The angle between the BLUE missile flight path and 
the target's direction of movement. 


CS 


The target course. 


DELTA 


The reciprocal of BETA. 


DELTAX 


Used for a) the X distance from the target ship to 
the RED missile, b) the present distance in the X 
direction from the firing ship to the target, and 
c) the distance in the X direction from the firing 
ship to the intercept point. 


DELTAY 


Same as DELTAX but in the Y direction. 


DIS 


Range to the intercept point. 
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DLT 


Range to the intercept point . 


DX 


D'istance in the X direction from the initial target 
position to the intercept point. Used if inter- 
cept point is over the firing ship. 


DY 


Distance in the Y direction from the initial target 
position to the intercept point. Used if the in- 
tercept point is directly over the firing ship. 


EX 


The exponent that defines the super ellipse of the 
intercept contour. 


FT 


BLUE missile flight time. 


FZL 


Left hand limit of the firing zone. 


FZL1 


Left hand limit of the second firing zone if a 
second exists. 


FZR 


Right hand limit of the firing zone. 


FZR1 


Right hand limit of the second firing zone. 


GAMMA 


Angle between the BLUE missile flight path and the 
line of sight to the target at launch. 


GSINE 


The sine of GAMMA. 


ID 


RED missile type. 


IND 


Flag indicating the feasibility of the intercept. 


IR 


RED missile's launching raid number. 


JOAN 


RED missile's target ship number. 


JULIA 


RED missile number for table lookup. 


K 


Firing ship number. 


KL 


Starting index for the BLUE missile class. 


KR 


Raider class number. 


LOOT I M 


Subroutine for determining raider position at a 
desired time. 


M 


Target number. 


NMS 


BLUE missile system number. 


NS 


The number of ships in the BLUE force. 
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PHI 


The BLUE firing bearing. 


PHI1 


The angle of axis rotation for intercept contour 
engagement evaluation. 


PI 


3.141592654. 


PIZ 


The hazard range squared. 


PK 


The single shot kill probability. 


PM 


The probability of successful launch. 


PROB 


PK times PM. 


RADIAN 


A function for converting angles in degrees to 
angles in radians. 


RAID 


The table with the raider position and weapon 
information . 


RAT 


A REAL*8 variable for use with RADIAN. 


REASON 


A code to indicate the reason for non-engagement. 


RMAX 


The maximum range of the BLUE missile. 


RMIN 


The minimum range of the BLUE missile. 


RUTH 


Target altitude in miles. 


RXATA 


Dummy vector use to determine the firing zone. 


SEEK 


Y intercept of the intercept contour seeker limi- 
tation approximation. 


SLOPE 


The slope of the line that approximates the seeker 
look limitation. 


SSPD 


Speed of the BLUE missile. 


THETA 


The angle, measured in the normal X,Y plane, of the 
line from the firing ship to the present target 
position . 


TH1 


The reciprocal of THETA. 


TIME 


The time of the proposed missile launch. 


TKIL 


The time of intercept. 


TX 


The time the target will enter the maximum missile 
range for this system. 
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TXL 


The time the target will depart the maximum missile 
range for this system. 


VEL 


The speed of the target. 


XI 


The X coordinate of the intercept. 


XIP 


The X intercept coordinate after axis rotation and 
translation . 


XIS 


The X intercept coordinate after axis translation. 


XP 


The X coordinate of the present target position. 


XT 


The end of game X coordinate for the target if it 
is a raider, or the target ship coordinate if the 
BLUE missile's target is a RED ASM/SSM. 


XW 


The initial position of the target in the X direc- 
tion . 


YI 


The Y coordinate of the intercept point. 


YIP 


The Y intercept coordinate after axis rotation and 
translation . 


YIS 


The Y intercept coordinate after axis translation. 


YP 


The Y coordinate of the present target position. 


YT 


The end of game Y coordinate for the target. 


YW 


The Y coordinate of the initial target position. 


ZLNKA 


A sorting subroutine. 



14. SUBROUTINE JAMM 



AREA 


Fourth root of target's radar cross section. 


COFF 


Jammer's main beam/side-lobe coefficient. 


DELX 


Target's track change in X direction minus 0.1 nm. 


DELY 


Target's track change in Y direction minus 0.1 nm. 


HJBW 


One-half jammer beam width in radians. 


HRBW 


One-half radar beam width in radians. 


ITS 


Target ship number. 


JJTBR 


Bearing from jammer to the jamming target. 
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JRBRG 


Bearing from the jammer to the detecting radar. 


LJBW 


Left bearing of the jamming beam. 


LTBRG 


Left bearing of the detecting radar's beam while 
pointed at the target. 


NASM 


Number of the ASM. 


NR 


Number of raiders. 


PJ 


Power of the jammer on the detecting radar's fre- 
quency . 


RCOF 


Detecting radar's side-lobe coefficient. 


RCOFF 


Detecting radar's side-lobe coefficient for specific 
jammer-target-radar geometry. 


RCPS 


Range to track CPA squared. 


RJBRG 


Detecting radar to jammer bearing. 


RJBW 


Right bearing of jammer beam. 


RJHS 


Detecting radar- jammer horizon squared. 


RRJS 


Range from t he detecting radar to the jammer 
squared. 


RTBRG 


Right bearing of the detecting r-adar's beam width 
while pointed at the target. 


SUMJP 


Sum of the jamming power arriving at the detecting 
radar . 


TBRG 


Target bearing from the detecting radar. 


TIMC 


Time constant used to adjust midpoint time. Base 
on target's launch time or game time of target's 
initial position. 


TIMI 


Time target is at its midpoint in detecting radar's 
detection envelope. 


TIMX 


Time target is at the X coordinate of its midpoint 
in the detecting radar's detection envelope. 


TIMY 


Time target is at the Y coordinate of its midpoint 
in the detecting radar's detection envelope. 


VD 


Divisor used in obtaining the X and Y components 
of the target's velocity. 
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vx 

VY 

XC 

XCOF 

XJ 

XJJ 

XJJT 

XJR 

XJT 

XO 

XR 

XRT 

XI 

YC 

YJ 

YJJ 

YJJT 

YJR 

YJT 

YO 

YR 

YRT 



X component of the target's velocity. 

Y component of the target's velocity. 

X coordinate of target's CPA to detecting radar. 

Product of the jammer and detecting radar coeffi- 
cients based on the jammer- target-radar geometry. 

X coordinate of the jammer when the target is at 
its midpoint. 

X coordinate of the jammer at game time 0. 

Difference in X coordinates of the jammer's posi- 
tion and the jammer's target position. 

Difference in X coordinates of the detecting radar's 
position and the jammer's position. 

X coordinate of the jammer's target. 

X coordinate of the target's terminal position. 

X coordinate of the detecting radar's position. 

Difference in Y coordinates of the detecting radar's 
position and the target's midpoint position. 

X coordinate of the target's initial position. 

Y coordinate of target's CPA to detecting radar. 

Y coordinate of the jammer when the target is at 
its midpoint. 

Y coordinate of the jammer at game time 0. 

Difference in Y coordinates of the jammer's posi- 
tion and the jammer's target position. 

Difference in Y coordinates of the detecting radar's 
position and the jammer's position. 

Y coordinate of the jammer's target. 

Y coordinate of the target's terminal position. 

Y coordinate of the detecting radar's position. 

Difference in Y coordinates of the detecting radar's 
position and the target's midpoint position. 
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Y1 Y coordinate of the target's initial position. 

Z1 Target altitude in nautical miles. 

15. SUBROUTINE BIDET 

A Coefficient used in solution of bistatic detection 

and radar horizon equations. See Appendix A for 
bistatic equation. 

A1 2A . 

A2 A squared. 

B Similar to A. See Appendix A. 

C Similar to A. See Appendix A. 

CC Constant used in solution of receiver and trans- 

mitter radar horizons equations. 

CHAR Constant of bistatic detection equation. See 

Appendix A. 

CHAS Bistatic radar detection range in clear environ- 

ment against one square meter target. 

COF Array used only in calling subroutine DTIMD. 

COFF See definition in subroutine JAMM portion of this 

glossary . 

CPRS Range squared from track CPA to bistatic receiver 

position . 

CPXS Range squared from track CPA to bistatic transmit- 

ter position. 

D Similar to A. See Appendix A. 

DELX Same definition as in subroutine CPA. 

DELY Same definition as in subroutine CPA. 

E Similar to A. See Appendix A. 

EE Constant used in the solution of bistatic receiver 

and transmitter radar horizon equations. 

HJBW Defined in subroutine JAMM. 

HORRS Receiver's radar horizon range squared. 
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HORXS Transmitter’s radar horizon range squared. 

HRBW Defined in subroutine JAMM. 

IFLAG Flag set to 1 for stationary targets. 

IT Ship number of bistatic transmitter. 

JAM Array indicating the number of active jammers on 

each radar frequency band. 

JFLAG Flag set to 1 when the algorithm has completed the 
jamming section of the detection phase. 

JJTBR Defined in subroutine JAMM. 

JRBRG Defined in subroutine JAMM. 

LJBW Defined in subroutine JAMM. 

LTBRG Defined in subroutine JAMM. 

NASM ASM number . 

PJ Defined in subroutine JAMM. 

Q Track change in Y direction divided by track change 

in X direction. 

RCOF Defined in subroutine JAMM. 

RCOFF Defined in subroutine JAMM. 

RD2 Square of the range from the bistatic receiver to 

stationary target's position. 

RJBRG Defined in subroutine JAMM. 

RJBW Defined in subroutine JAMM. 

RJHS Defined in subroutine JAMM. 

RRAD Radical part of solution to quadratic equation for 

time above bistatic receiver's radar horizon. 

RRH Range to receiver's radar horizon. 

RRJS Defined in subroutine JAMM. 

RTBRG Defined in subroutine JAMM. 

SUMJP Defined in subroutine JAMM. 
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T 

TBRG 

TIMC 

TIME 

TIMEI 

TIMI 

TR1 

TR2 

TX1 

TX2 

T1 

T10 

T2 

T20 

V 

VD 

VX 

VY 

XCOF 

XCOFF 

XD2 

XJ 

XJJ 



An array containing the bistatic detection times. 
Defined in subroutine JAMM. 

Time constant based on target's launch time of 
missile. Otherwise start of game time. 

Array containing the real parts to the solution of 
the bistatic equation. 

Array containing the complex parts to the solution 
of the bistatic equation. 

Defined in subroutine JAMM. 

Time target enters receiver's radar horizon range. 

Time target departs receiver's radar horizon range. 

Time target enters transmitter's radar horizon 
range . 

Time target departs transmitter's radar horizon 
range . 

Time target enters bistatic radar detection envelope. 

Time target is first above both transmitter and 
receiver's radar horizon. 

Time target departs bistatic radar detection en- 
velope. 

Time target is no longer above both the transmitter 
and receiver's radar horizon. 

Target's velocity in nautical miles per minute. 
Defined in subroutine JAMM. 

X component of target velocity. 

Y component of target velocity. 

Defined in subroutine JAMM. 

Defined in subroutine JAMM. 

Square of the range from the bistatic transmitter 
to the stationary target's position. 

Defined in subroutine JAMM. 

Defined in subroutine JAMM. 
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XJJT 



Defined in subroutine JAMM. 



XJR 

XO 

XR 

XRAD 

XRC 

XRT 

XX 

XXC 

XI 

YJ 

YJJ 

YJJT 

YJR 

YJT 

YO 

YR 

YRC 

YRT 

YX 

YXC 

Y1 

ZR 

ZRN 

ZX 



Defined in subroutine JAMM. 

X coordinate of target's terminal position. 

X coordinate of bistatic radar receiver position. 
Similar to RRAD above. 

X coordinate of target track CPA to bistatic 
receiver . 

Defined in subroutine JAMM. 

X coordinate of the bistatic transmitter. 

X coordinate of target track CPA to bistatic 
transmitter . 

X coordinate of target's initial position. 

Defined in subroutine JAMM. 

Defined in subroutine JAMM. 

Defined in subroutine JAMM. 

Defined in subroutine JAMM. 

Defined in subroutine JAMM. 

Y coordinate of target's terminal position. 

Y coordinate of bistatic receiver's position. 

Y coordinate of target track CPA to bistatic 
receiver . 

Defined in subroutine JAMM. 

Y coordinate of the bistatic transmitter position. 

Y coordinate of target track CPA to bistatic trans 
mitt er . 

Y coordinate of target's initial position. 

Bistatic receiver antenna height in feet. 

Bistatic receiver antenna height in nautical miles 
Bistatic transmitter antenna height in feet. 
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ZXN Bistatic transmitter antenna height in nautical 

miles . 

Z1 Target altitude in nautical miles. 

Z12 Target altitude in nautical miles squared. 



16. SUBROUTINE CPA 



DELX 



DELY 

Q 

XC 

XI 

XO 

XP 

YC 

YI 

YO 

YP 



First used as a test variable to determine if track 
passes through reference position. Second usage is 
to indicate track change in X direction. 



First used as a test variable to determine if track 
passes through reference position. Second usage is 
to indicate track change in Y direction. 

Track change in Y direction divided by track change 
in X direction. 



X coordinate of 
X coordinate of 
X coordinate of 
X coordinate of 

Y coordinate of 

Y coordinate of 

Y coordinate of 

Y coordinate of 



CPA position, 
initial track position, 
terminal track position, 
reference position. 

CPA position, 
initial track position, 
terminal track position, 
reference position. 



17. SUBROUTINE RCHECK 

DELX The distance along the X axis from the raider 

position to the potential target position. 

DELY The distance along the Y axis from the raider 

position to the potential target position. 

DIS The square of the distance from the raider to the 

potential target. 

1IT The height of the highest radar on the potential 

target ship. 

IITR The radai’ horizon for IIT. 



180 



ID 


The class number of the potential target ship. 


IR 


The raider number. 


IT 


The potential target ship number. 


K 


Flag for go/nogo situation. 


RADAR 


The combined radar horizons of the raider and the 
potential target. 


RADAR2 


RADAR squared . 


RR 


The maximum theoretical range for the raider's 
radar . 


RR2 


RR squared. 


RTR 


Radar horizon for the raider. 


SHIP 


The table containing the BLUE ship location infor- 
mation. 


XS 


The raider's present position along the X axis. 


XT 


The ship's present position along the X axis. 


YS 


The raider's present position along the Y axis. 


YT 


The ship's present position along the Y axis. 



18. SUBROUTINE TMDASM 



A 


Coefficient A in the quadratic equation. 


B 


Coefficient B in the quadratic equation. 


B1 


B in Y = mX + b 


C 


Coefficient C in the quadratic equation. 


DELTAX 


The distance in the X direction from present posi- 
tion to the end of game position. 


DELTAY 


The distance in the Y direction from present posi- 
tion to the end of game position. 


DELX 


The distance in the X direction from the intercept 
position to the present position. 


DELY 


The distance in the Y direction from the intercept 
position to the present position. 
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DEX 


The distance in the X direction from the firing 
ship to the target at present time. 


DEY 


The distance in the Y direction from the firing 
ship to the target at present time. 


DIC 


The distance in either the X or Y direction from 
intercept to end of game. 


DIF 


The distance in either the X or Y direction from 
present position to end of game position. 


Dl 


Distance from present position to first intercept 
point, squared. 


D2 


Distance from present position to second inter- 
cept point, squared. 


IFLAG 


Flag indicating target is already within the indi- 
cated range. 


IRAID 


Raider number of target or launching raid. 


ISAMBT 


Battery number. 


ITS 


Target ship number when target is an ASM/SSM. 


KFLAG 


Flag indicating target moving in a north or south 
direction . 


M 


Raider number. 


NFS 


Firing ship number. 


RAID 


The raid array. 


RIP 


The range to the target or to intercept, whichever 
is smaller. 


ROOT 


B 2 - 4AC. 


SHIP 


The ship array. 


SLM 


m in Y = mX + b. 


SRG 


Range of interest. 


TIME 


Time target crosses range inbound. 


TIML 


Time target crosses range outbound. 


TML 


Present game time. 


TOAST 


B 2 - 4AC . 
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VEL 


Target speed. 


XI 


Inbound crossing point. 


YI 


Inbound crossing point. 



19. SUBROUTINE TMDBOM 



A 


A in the quadratic equation. 


B 


B in the quadratic equation. 


BA 


b in Y = mX + b . 


C 


C in the quadratic equation. 


DELTAX 


The distance from present position to the end of 
game position in the X direction. 


DELTAY 


The distance from the present position to the end 
of game position in the Y direction. 


DELX 


Distance from the intercept point to present posi 
tion in the X direction. 


DELY 


Distance in the Y direction from the intercept 
point to the present position. 


DIC 


Distance from the intercept point to the end of 
game position in the X direction. 


DIF 


Distance from the present position to the end of 
game position in the X direction. 


IR 


Raider number. 


RAID 


Raid array. 


RGE 


The launch range. 


SLM 


m in Y = mX + b. 


T 


Time raider will reach the launch range. 


VEL 


The raider speed. 
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20. SUBROUTINE LOCTIM 



D The distance between the raider end of game posi- 

tion and the position stored in the RAID table. 

DELTAX The distance along the X axis from the end of game 
position to the position stored in the RAID table. 

DELTAY Distance along the Y axis. 

DL Distance traveled in the time interval defined as 

the difference between the time stored in the RAID 
table and the time for which the raider position i 
desired . 

K The raider number. 



PTIME The time stored in the RAID table. 

RAID The table that contains the position and weapon 

information . 

T The time for which the raider position is desired. 

VEL The raider speed. 

XP The X position at time T. 

XZ The X position stored in the RAID table. 

YP The Y position at time T. 

YZ The Y position stored in the RAID table. 



21. SUBROUTINE ASMTIM 

ASMT The table containing the data on the RED missile 

engagement . 

D The distance between the missile launch point and 

the impact point. 

DL The distance the missile will travel in the time 

defined as the difference between the launch time 
and the time for which the missile position is 
desired . 

ID The RED missile type. 

M The RED missile target number. 

MA The RED missile number for table lookup. 
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T 

VEL 

XP 

YP 



The time for which the ASM/SSM position is desired. 
The missile speed. 

The X position at the time T. 

The Y position at the time T. 
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210 WRITE (6,570) K JLE 1720 

NADAR(I,K) = 0 JLE 1730 

220 CONTINUE JLE 17^0 

JLE 1750 

230 IF (KF0G.EQ.1) GO TO 270 JLE 1760 
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